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I. Introduction 

The great importance of impurity states in determining the electrical 
properties of semiconductors has been recognized for a long time.^ There 
are two varieties of these states, the so-called donor and acceptor states. 
In a donor state an electron is localized around an impurity from which 




Empty conduction 
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Donor levels 


Acceptor levels 
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Fig. 1. Donor and acceptor levels in the band gap of a semiconductor (schematic). 



it can be released into the conduction band with a rather small expendi- 
ture of energy. In this way the impurity atom (donor) can make an elec- 
tron available for electrical conduction (Fig. 1). Similarly, an acceptor 
state is an available electronic state around an impurity which can accept 
an electron from the valence band. When this transition occurs an elec- 
tron is taken from the valence band, or a hole is produced (Fig. 1). An 
electrical current can now be carried by the valence band. Since the en- 
ergies involved in these processes are of the order of 0.01 to 0.1 ev, and 
hence are in general much smaller than the energy required to excite an 
electron directly from the valence band into the conduction band (of the 
order of 1 ev), the electrons and holes which are thermally excited from 
impurity states generally dominate the electrical properties of semicon- 


1 N. F. Mott and R. W. Gurney, Electronic Procession in Ionic Crystals,” p. 166. 
Oxford Univ. Press, London and New York, 1940. G. L. Pearson and J. Bardeen, 
Phya. Rev. 75, 865 (1948). Previous reviews dealing with impurities in silicon and 
germanium are: J. A. Burton, Phyaica 20, 845 (1954) and E. Burstein, G. S. Picus, 
and N. Sclar, m “Photoconductivity Conference held at Atlantic City, Nov. 4-6” (R. 
G, Breckenridge, ed.), p. 353, Wiley, New York, 1954. 
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ductors at all but the highest temperatures. Quite apart from their role as 
sources of charge carriers, however, impurity states have become the ob- 
ject of both experimental and theoretical studies in their own right, and 
their properties have been mapped out and interpreted in considerable 
detail in recent years. 

Most progress in this field has been accomplished by studying ger- 
manium and silicon whose impurity content can now be controlled over 
an impressive range extending from about one part in 10® to one part in 
10®. The high purity attainable in these materials has also made it pos- 
sible to elucidate the electronic structure of the host crystals in an espe- 
cially detailed way. This in turn has enabled one to interpret the proper- 
ties of the impurity states more fully. For this reason we shall, in this 
article, focus our attention principally on impurity states in silicon and 
germanium. 

Some impurities in silicon and germanium, such as Cu and Au, can 
bind carriers with quite substantial energies, of the order of 0.5 ev. The 
orbits associated with these states are presumably rather small, of the 
order of a few angstroms and their structure is not very well understood. 
Apart from giving a brief summary of their properties in the Appendix, 
we shall not discuss such ‘^deep” states further in the present article. 
Instead we shall limit ourselves to the ‘^shallow” states, whose orbits 
are of the order of 50 A in diameter, and whose binding energies are be- 
tween about 0.01 to 0.1 ev. 

It may seem very surprising at first sight that the orbit of a trapped 
electron or hole wending its way through hundreds of crystal cells can be 
described at all simply. The theory shows, however, that the larger the 
orbit the more accurately it can be understood! 

To appreciate this point let us consider a donor state for a moment. 
Such a state can be pictured roughly as follows. Assume a neutral impurity 
atom is introduced into the lattice and imagine that one electron is re- 
moved from it, leaving a positive ion behind. This ion will polarize the 
semiconductor so that at large distances it produces an electrostatic 
potential 

= — ( 1 ) 

Kr 

Here k is the static dielectric constant. When the electron is brought 
back, one of two things may happen. It may be energetically favorable 
for this electron to occupy an orbit rather near the impurity ion. In this 
case Eq. (1) does not apply and we have a rather complicated state of 
affairs in the immedi^e vicinity of the impurity ion. In brief, we obtain 
a ^Meep^’ impurity state. On the other hand, the lowest energy state 
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may have a large orbit ov^r most of which Eq. (1) appUes. One would 
then guess that such a state can be described by the Schroedinger equation 

where m* is an appropriate effective mass,^ 

It is clear that Eq. (2) will apply more accurately the larger the orbit, 
since then the region near the impurity ion, where Eq. (1) does not hold, 
plays an insignificant role. Here we see one of the main reasons why we 
can form an accurate picture of the shallow impurity states: The poten- 
tial is basically Coulomb-like and hence the impurity states are very 
closely related to the states of the hydrogen atom. A second reason is 
slightly more subtle. We shall see later that the concept of an effective 
mass holds only for a state whose wave function F(r) has a slow spatial 
variation. This principle also favors the description of shallow impurity 
states with large orbits. 

The analogy with the hydrogen atom suggests the existence not only 
of a lowest bound state, but of a whole spectrum of excited bound states. 
Such states have been clearly identified experimentally.® Since their 
orbits are even larger than that of the ground state, equations of the type 
(2) are particularly applicable to them. 

Our understanding of the physical properties of impurity states is 
still being rapidly expanded. However, a great deal of insight into these 
fascinating “ atoms has already been gained. They are very delicate 
structures. They may be disrupted by an energy of the order of 0.01 ev 
and their orbits have the relatively immense dimensions in the range from 
10”^ to 10“* cm. Nevertheless, their structure is determined by the same 
quantum-mechanical laws and the same Coulomb attraction which gov- 
ern their robust counterparts, the mesic atoms, with binding energies 
of about 10* ev and orbital diameters of about 10”^^ cm. The range of 
energies, dimensions, and physical conditions in which hydrogen-like 
atoms occur in nature is truly remarkable. 

The plan of the present article is as follows. First we shall review the 
available experimental information concerning the shallow impurity 
states in silicon and germanium (Part II, Sections 1-3). This information 
has been obtained by studying the electrical properties of sihcon and 
germanium containing impurities, by infrared absorption experiments, 
by electron and nuclear spin resonance experiments, and by magnetic 
susceptibility measurements. Following this, we shall present the theory 

* The actual band structures of silicon and germanium are rather complex, which 
requires a modification of the kinetic energy term (see Sections 5 and 9). 

® E. Burstein, E. E. Bell, J. W. Davisson, and M. Lax, Phys. Chem. 67, 849 (1953), 
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which has been developed to account for the empirical facts and shall 
discuss its successes and failures (Parts III-VI, Sections 4~15). Finally 
the properties of the deep-lying impurity states in silicon and germanium 
will be summarized in the Appendix. Throughout this article we shall 
discuss only isolated impurity states for which overlap between nearest 
neighbors can be neglected. 

II. Empirical Properties 

In this part we shall describe briefly the most important empirical 
information concerning the shallow impurity states in silicon and ger- 
manium. We shall return to these data in subsequent parts where they 
will be discussed in the light of the theoretical developments. 

1. Energy Levels 

a. Ionization Energies 

The ionization energy of a shallow impurity ^^atom^’ is the energy of 
its ground state relative to that of the nearest band edge; or, physically 
speaking, the minimum energy required to liberate a bound carrier from 
its impurity ion. We shall denote this (positive) energy by 

The ionization energy has been determined empirically in two ways. 
In one of these the number n of thermally released, conducting carriers 
is measured by means of the Hall constant, which is inversely proportional 
to n. Since the temperature dependence of n is governed primarily by an 
appropriate Boltzmann factor, which under the usual experimental con- 
ditions is exp {—Eo/kT)y^ this is a convenient method of determining 
Eo- Early measurements of this kind were made by Debye and Conwell^ 
and an extensive analysis of Hall effect data has been carried out by the 
Bell Telephone Laboratories group.®-^*® The results are listed in Table I 
under the heading ‘‘Thermal.^’ 

The other method makes use of infrared irradiation at low tempera- 
tures. One observes the photon energy required to raise a carrier from the 
ground state to an excited bound state and adds to this the theoretically 
calculated energy of the excited state. It is believed that the latter, which 
is a small fraction of the total ionization energy, can be calculated quite 

* W. Shockley, “Electrons and Holes in Semiconductors,” Section 16.4. Van Nostrand, 
New York, 1950. 

* P. P. Debye and E. M. Conwell, Phys. Rev. 98 , 693 (1954). 

* T. H. Geballe and F. J. Morin, Phys. Rev. 96, 1085 (1954). Due to an oversight, the 
temperatures in their graph are too high by a factor of 10. 

' F. J. Morin, J. P. Maita, R, G. Shulman, and N. B. Hannay, Phys, Rev. 96 , 833 
(1953). The value of 0.039 ev given there for P is a misprint. 

* F. J. Morin, private communication. 
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reliably. Such studies have been carried out particularly by Burstein and 
his co-workers. Their results are given in Table I under the heading 
‘'Optical.'^ An alternative infrared technique is to determine the thresh- 
old energy for photoconductivity. Such experiments have been per- 
formed® but yield only very rough estimates of Eq, 

Table I. Ionization Energies 


(ev)“ 


Impurity element Thermal Optical 


0.045(11) 
0.053 (10) 
0.043 (10) 

0.046 (9) 
0.067 (9) 
0.071 (9) 
0.154 (9) 



p ] 


0.0120(6) 


As 

• Donors 

0.0127 (6) 


SbJ 


0.0096 (6) 

In germanium 


I 

0.0104(6) 


A1 1 
Gai 

> Acceptors 

0.0102 (6) 

0.0108 (6) 


In J 

1 

0.0112 (6) 


In silicon 



0.033 (7) 
0.044 (7) 
0.049 (7) 
0.039 (7) 
0.069 (8) 
0.045 (7) 
0.057 (7) 
0.065 (7) 
0.16 (7) 


« Numbers in parentheses refer to footnote references in text. 

It will be seen from Table I that, even though the results of the ther- 
mal and optical methods are in fair agreement, there remain discrepan- 
cies which, in some cases, are 10 to 20%. It appears very unlikely that the 
deviations are due entirely to experimental error. 

The following trends are apparent from this table. 

1. All shallow donors (except Li) are group V elements; all shallow 
acceptors are group III elements.^®® 

* E. Burstein, G. S. Picus, B. Henvis, and R. Wallis, /. Phys. Chem. Solids 1, 65 (1956). 

G. S. Picus, E. Burstein, and B. Henvis, J. Phys. Chem, Solids 1, 75 (1956). 
loo suggests the following “chemical*^ models: A group V donor enters the lattice 
substitutionally. It contributes four of its five valence electrons to saturate all 
covalent bonds. The fifth electron is weakly bound by the Coulomb field of the 
resulting positive donor ion. Similarly, a group III acceptor enters the lattice sub- 
stitutionally. Its three valence electrons saturate all but one of the covalent bonds. 
The remaining ^^hole’’ in the bond structure circulates in a large orbit around the 
net negative charge in the vicinity of the acceptor ion. The energy required to 
remove this hole to infinity is the ionization energy of the acceptor. 

“ G. S. Picus, private communication. These data supersede those quoted in ref. 10. 



Table II. Spectra of Excited States in Silicon 


Impurity element 


Energies relative to ground state® 
(in units of 0.01 ev) 


p 


3.45 (11) 

3.95 (10) 

4.26 (11) 

4.46(10) 



As 

Donors 

4.21 (10) 

4.74 (10) 

5.06 (10) 

5.29(10) 



Sb) 


3.18 (11) 

3.65 (10) 

3.99 (10) 




B 


3.0 (13) 

3.5 (13) 

3.85 (13) 

3 97 (13) 

4.14 (13) 

4.21(13) 4.27(13) 

AJ 1 

Acceptors 

5.50 (13) 

5.87 (13) 

6.43(13) 

6.51 (13) 

6.72 (13) 

6.84 (13) 

Gal 

5.85(13) 

6.22 (13) 

6.43 (13) 

'-O.S (13) 

6.95(13) 

7.25 (13) 

In J 


14.1 (10) 

14.6 (12) 

15.1 (12) 





® Numbers in parentheses refer to footnote references in text. 
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2. The ionization energies of all donors and acceptors are very similar 
in germanium, about 0.011 ev. 

3. The ionization energies in silicon are larger than in germanium 
and show a much wider spread, varying from 0.03 to 0.16 ev. 

The ionization energies are discussed further in Sections 6, 7, and 10. 

6. Spectra of Excited States 

Transitions of carriers from the ground state of the impurity to dis- 
crete excited levels were first observed by Burstein and coworkers^ with 



Fig. 2. Infrared absorption lines of B, Al, and Ga acceptors in silicon. T is the 
intensity of the transmitted radiation. (Courtesy H. J. Hrostowski,^*) 

the use of infrared absorption. Further studies of such transitions have 
been carried out since by R. Newman, the group at the Naval Re- 
search Laboratories*’^® and H. Hrostowski and R. H. Kaiser.^® The 
excited states are of great theoretical interest. As a consequence of their 

R. Newman, Phys, Re&, 99 , 465 (1955). 

R. Newman, Pkys. Reo. 108, 103 (1956). 

H. J. Hrostowski and R. H. Kaiser, BulL Am. Phys. Soc. [2] 2, 66 (1957); also 
H. J. Hrostowski, private communication. 
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very large orbits, they are insensitive to the complex state of affairs in 
the immediate vicinity of the impurity ion and should, therefore, be 
described well by an effective mass theory. Table II gives the experimen- 
tal values of the energies of the absorbed photons, that is, the differences 
in energy between the ground state and the various excited states. Only 
impurity states in silicon have been investigated by this method so far 
because the radiation needed in the case of germanium lies in an experi- 
mentally difficult region (X > 100 m)* Figure 2 shows typical absorption 


CONDUCTION 



Fig. 3. Energy levels of donor states in silicon, experiment and theory. The posi- 
tion of the conduction band edge relative to the bound levels cannot be fixed very 
accurately by experiment and is here taken at 0.0109 ev above the first excited state, 
as predicted by theory. 

plots, obtained by Hrostowski.^® The positions of the excited levels rela- 
tive to each other are of greater significance than the positions of the 
excited levels relative to the ground state. They are displayed in Figs. 
3 and 4. The figures also show the positions of the ground states and, 
approximately, those of the nearest band edge. We note a striking simi- 
larity between the excited levels of different donors and, to a lesser ex- 
tent, between the excited levels of different acceptors. 

At helium temperature most of the lines have a width between half 
power points in the range from 0.5 to 1.0 X 10”^ ev.^® At nitrogen tem- 
perature the width increases to about 2 X 10“^ ev in the case of On 
the other hand, the lines associated with In broaden almost to the point of 
M. Lax and E. Burstein, Phys. Rev, 100, 592 (1955). 
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band edge is here taken to be 0.014 ev below the first excited state. The resulting 
ground state energies differ slightly from those quoted in Table I. 

disappearing^^* at nitrogen temperature. The line widths are interesting 
because they reflect the interaction of the bound carriers with the lattice 
vibrations. They are discussed further in Section 15. 

2. Spin Resonance 

a. Electron Spin Resonance 

The spin resonance of bound donor electrons in silicon was first ob- 
served by Fletcher et and has since been extensively studied.^^*^® 

This behavior is probably not typical of shallow impurity states and connected 
with the large ionization energy (0.16 ev) of In acceptors. 

R. C. Fletcher, W. A. Yager, G. L. Pearson, A. N. Holden, W. T. Read, and F. R. 
Merritt, Phyi. Rev. 94, 1392 (1954), 

R. C. Fletcher, W. A. Yager, G. L. Pearson, and F, R. Merritt, Phya. Rev. 96 , 844 
(1954). 

A. Honig and A. F. Kip, Phya. Rev. 96 , 1686 (1954); A. Honig, Phya. Rev. 96 , 234 
(1954). 

G. Feher, Phya. Rev. 103, 834 (1956) and private communication. 
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(No resonance has yet been observed either for bound electrons in ger- 
manium or for bound holes in silicon and germanium.) 

These measurements have thrown a great deal of light on the wave 
functions of the donor states, for, as we shall see, the resonance frequen- 
cies, the line widths, and the spin-lattice relaxation times all depend 
crucially on the hyperfine interaction of the donor electron with the 
nuclear magnetic moments within its orbit. The hyperfine interaction, 
in turn, gives us information about the square of the electronic wave 
function at the positions of the nuclei. 

In their original experiments, Fletcher and co-workers found that the 
resonance pattern of silicon containing P donors consisted of an extremely 
well-resolved doublet (width « 3 oersted, separation « 42 oersted) ; the 
pattern of As was a quartet; and that of Sb was a superposition of a 
strong sextet and a weaker octet. When these observations were coupled 
with the fact that the naturally occurring nuclei of P have spin those 
of As have spin f, and those of Sb have spins i and it was clear that the 
observed structures were the result of the hyperfine interaction of the 
donor electrons with their donor nuclei. At fixed frequency, the separa- 
tion of successive lines arising from this interaction is given by 

A// =yJ^lf(0)P, (2.1) 

where hd and I jy are the magnetic moment and spin of the donor nucleus, 
respectively, and |\^^(0)|^ is the probability density of the donor elec- 
tron at the position of the donor nucleus. Thus the measurements of 

Table III. Probability Densities at Donor Nuclei in Silicon'* 


Donor 

|^(0)i* (cm'*) 

Li 

0.0032 X 10*M18) 

P 

0.44 X 10«<16) 

As 

1.80 X 102^(16) 

Sb 

1.20 X 10*" (16) 


» Numbers in parentheses refer to footnote references in text. 

AH yield the important quantity i^(0)12. The results are given in Table 
III. The table also includes the value of |^(0)1^ for Li obtained by Feher 
using a double resonance technique (see Section 2b). A trace of the elec- 
tron resonance lines obtained with P donors is shown in Fig. 5. 

Each resonance line of the hyperfine multiplets was found to have a 
width between half power points of about 3 oersteds. Fletcher et 
observed that this width could be associated with so-called inhomoge- 
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neous broadening^® which suggested strongly that it originated in further, 
unresolved hyperfine interactions of the electron. Now 4.7 % of natural 
silicon is present in the form of the isotope Si^®, which carries a magnetic 
moment. Therefore, the resonance field of the electron should depend not 



Fig. 5. Electron spin resonance lines of P donor states, at a frequency of 9000 me. 
The doublet corresponds to the two possible orientations of the P nuclear spin. The 
width of each line is the result of the interactions of the electron with the moments 
of Si** nuclei. (Courtesy G. Feher.**) 

only on the orientation of the spin of the donor nucleus, but also on the 
position and orientation of each Si*® spin within its orbit. Since the inter- 
action with the Si*® moments is rather weak — as we shall see in more 
detail in Section 7c, — it offered a natural explanation of the observed 
widths. Table IV contains a list of the empirical line widths of the spin 
resonances of various donors. 

Table IV. Widths of Donor Resonances Due to Si” 


Donor 

Width between half-value points 
of absorption line (oersted)* 

Li 

3.0(20) 

P 

2.9(20) 

As 

3.6 (20) 

Sb 

2.7 (20) 


“ Numbers in parentheses refer to footnote references in text. 

Definitive experimental results concerning the spin lattice relaxa- 
tion times of isolated donor atoms, are still lacking. The observed times 
are very long at low temperatures and, therefore, extremely sensitive to 
small disturbances such as neighboring impurities or a small number of 
conduction electrons. There are two types of relaxation processes:*^ 1. 
Those involving the simultaneous flip of an electron spin and a nuclear 

” A. M. Portis, Phys. Rev. 91, 1071 (1955). 

G. Feher, private communication. 

D. Pines, J. Bardeen, and C. P. Slichter, Phys. Rev. 106, 489 (1957). 
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spin (donor or which are characterized by the time T*, and 2. 

those in which only the electron spin flips, which are characterized by the 
time T,. T, presumably originates in the hyperfine interaction, modu- 
lated by the lattice vibrations. The origin of T* is somewhat uncertain. 
To give some idea of the times involved, we quote results from the experi- 
ments of Abragam and Cambrisson.^^ Working with about 10^^ As 
atoms/cc, at 9000 me and 2^K, they find T® = 10 min and T, = 20 min. 

b. Double Resonance 

A new double resonance technique has been applied successfully by 
Feher^* to donor states in silicon. To illustrate the method we consider 

a' 



Fig. 6. Magnetic level scheme for a double resonance experiment on a P donor 
state (schematic) ; hv^ corresponds to an electron spin flip, hvn^ and hvn~ to nuclear 
spin flips. (After G. Feher.^*) 

the electron associated with a P donor (nuclear spin, / = i) and ignore 
for the moment the effect of the Si^® nuclei. When the system consisting of 
electron spin and nuclear spin is placed in a strong static magnetic field, 
it has the energy level spectrum shown in Fig. 6. The transition from A 
to A', corresponding to a microwave frequency, is associated with one 
line of the hyperfine doublet of the electron spin resonance; that from 
B to B' is associated with the other line. The separations hvn'^ and hvrTy 
which lie in the radio-frequency range, and correspond to nuclear spin 
flips are given by 

hPn^ = 2up ± y !*)• (2.2) 

** A. Abragam and J. Combrisson, Compt. rend. 243, 576 (1956). 
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Here fip is the magnetic moment of the P nucleus, H the external field, 
and ± (8ir/3)Mol^(0)|^ is the effective field at the nucleus arising from the 
electron spin. Suppose that the electron resonance hve is partially satu- 
rated. The signal arising from these transitions is proportional to the 
difference in population between A' and A, If the nuclear transition 
hvn’^ (or hvtT) is also induced the population difference is altered. Hence 
the signal strength of the electron spin resonance is changed. In this way 
the nuclear frequencies Vn'^ and VrT can be determined by observing the 
electron resonance as a function of the nuclear frequency. The value of 
|^(0)12 can then be determined using Eq. (2.2). 



Fig. 7. Double resonance. Electron spin resonance absorption as a function of the 
“nuclear” frequency vn. Each line corresponds to the resonance frequency of a Si** 
nucleus at a particular site relative to the P donor. (Courtesy G. Feher.*®) 


Feher was able by this method to determine |^(0)|^ for Li donors, 
which has the rather small value given in Table III. Ordinary techniques 
of electron resonance were not adequate for resolving the hyperfine 
multiplet. 

Another interesting application of the technique of double resonance 
involves the Si^® nuclei. Again let us take a P donor and consider a partic- 
ular Si®® nucleus at site within the orbit of the donor electron. Clearly 
the nucleus has resonance frequencies given by 

= 2ms, ± y (2.2') 

in a magnetic field, depending on the orientation of the electron spin; 
fisi is the magnetic moment of the Si®® nucleus. As before, a nuclear signal 
given by (2.2') changes the intensity of the electron absorption line if the 
electron spin resonance is partly saturated. In this way the values of 
|^(rz)|® can be obtained for all lattice sites Ti. Figure 7 shows the electron 
resonance signal as a function of the nuclear frequency v for the case of 
a P donor. Each line yields the value of at some lattice site rz. Thus in 
a sense the entire wave function of the donor electron is mapped. The 
only remaining problem is the assignment of the various observed fre- 
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quencies to specific lattice sites. We shall return to a more detailed anal- 
ysis of such experiments in Section 7c. 


3. Static Magnetic Susceptibility 

Impurity atoms exhibit both a temperature independent diamagnet- 
ism and a paramagnetism varying as l/T. As a very rough estimate we 
expect, per impurity atom, 


Xdia 


and 


. 2 \ f “ 10~^® cgs units in silicon .. 

6m*c^ ^ I — cgs units in germanium * 

0 23 

Xpara ^ Xfrae ap.o X CgS UOitS; (3.2) 


here m* is an average effective mass and T is measured in ®K. For com- 
parison, the diamagnetic susceptibilities of the host materials are 

Xsi = ^2.7 X 10“^ cgs volume units 

XGe = —5.6 X 10“^ cgs volume units. ^ 

Present techniques for measurement of the magnetic susceptibility 
have an accuracy of 1 or 2%. Therefore, to distinguish clearly the con- 
tributions of the impurity states requires concentrations of the order of 
10^® to 10^^ impurities per cm®. Unfortunately, overlap effects begin to 
set in just in this range. As a result, no firm conclusions about the mag- 
netic characteristics of isolated impurity states can as yet be drawn even 
though both diamagnetic and temperature dependent paramagnetic 
effects of the correct order of magnitude have been clearly established by 
experiment.^® It is to be hoped that techniques will soon be improved by 
a factor of about 3 to permit experimentation on more dilute samples. 
The theory of the magnetic susceptibilities is discussed in Section 13. 

III. Structure of Donor States 

4. Conduction Bands of Silicon and Germanium 

A donor state may be regarded as a wave packet which, as we shall see 
in the next section, is mainly composed of Bloch waves chosen from near 
the bottom of the conduction band. Thus a knowledge of the behavior 
of the conduction band near its lowest point (or points) is required for 
a quantitative understanding of these states. Fortunately, we have ac- 
quired a very detailed picture of the band structure of silicon and ger- 
manium, during the last few years, primarily thanks to the successful 

*3 A. Van Itterbeck and W. Duchateau, Physica 22, 649 (1956); also R. Bowers, 
private communication. 
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execution of cyclotron resonance experiments in these materials.^^ 
We shall summarize the relevant information obtained from these 
experiments. 

a. Silicon 

Figure 8 shows schematically the energy band structure of silicon 
along the (1,0,0) direction in k space.^® The point A represents one of the 
6 equivalent minima of the conduction band which are located at the 



Fig. 8. Energy band structure of silicon in the (1,0,0) direction (schematic). 

points (fco,0,0), ( — A:o,0,0), . . . (0,0,— fco). The exact magnitude of Ao is 
not yet known, but there are indications that the conduction band minima 
are about three-quarters of the way between k = 0 and the zone bound- 
aries.^® Two Bloch waves, corresponding to the two possible orientations 
of the electron spin, belong to each point such as A, Near one of these 
minima, say the one on the (0,0,1) axis, the energy of the Bloch waves 
is given by an expression of the form 


^ 2mj 2m, 

Here the zero of energy is taken at the bottom of the conduction band; 

For a general discussion see G. Dresselhaus, A. F. Kip, and C. Kittel, Phya, Rev. 98 , 
368 (1955). 

F. Herman, Phys. Rev. 96, 847 (1955); Physica 20, 801 (1954). 

G. G. McFarlane and V. Roberts, Phys. Rev. 98 , 1965 (1955); W. Kohn, Phya. Rev. 
98 , 1561 (1955) and unpublished work. 
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the longitudinal and transverse ‘‘masses” mi and have the following 
values 

mi — (0.98 ± 0.02)wo on 

mt — (0.19 ± 0.01)mo. ^ ' 

b. Germanium 

The energy band structure of germanium along the (1,1,1) direction 
in k space is shown schematically in Fig. 9.^® It is known with certainty 
from the behavior of the cyclotron resonance pattern that the conduction 
band minima lie on the (1,1,1) axis and the equivalent axes. There is 



Fig. 9. Energy band structure of germanium in the (1,1,1) direction (schematic). 

other strong evidence that the minima occur at the zone boundary,^* as 
shown in Fig. 9. Since points such as (1,1,1) and ( — 1,-1, — 1) are then 
equivalent, because they differ by a reciprocal lattice vector, there are 
four equivalent conduction band minima. If we take the (1,1,1) direction 
as a new z axis and denote the distance in k space of the point A from 
k = 0 by fco, the energy in the vicinity of A is again given by an expression 
of the form (4.1). The constants mi and mt have the following values in 
germanium:^® 

mi — (1.60 ± 0.008)mo .. on 

mt = (0.0813 ± 0.002)mo. ^ 

R. N. Dexter, B. Lax, A, F. Kip, and G. Dresselhaus, Phys, Rev. 96, 1222 (1954) 
and private communication from Dr. Zeiger. 

** D. K. Stevens, J. W. Cleland, J. H. Crawford, Jr., and H. C. Schweinler, Phys. Rev. 
100, 1084 (1955). 

R. C. Fletcher, W. A. Yager, and F. R. Merritt, Phys. Rev. 100, 747 (1955). 
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5. Effective Mass Theory of Donor States 

In this section we shall show that, with certain reservations, the wave 
function of a donor state in silicon or germanium can be expressed in the 
form 

N 

^(r) = 2 c ^ FAr )< pAr ). (5.1) 

1 


Here N is the number of equivalent minima, which is six in silicon and 
probably four in germanium; the aj are certain numerical coefficients; 
<Pj(T) is the Bloch wave at the 7th minimum; and the functions Fy(r) are 
^‘hydrogen-like^' envelope functions. If the energy of the conduction 
band near the jib. minimum has the form 


3 

E,ik) = ^ Di‘‘%akp (5.2) 

k being reckoned from the position of the mimimum, the functions Fy(r) 
satisfy so-called effective mass equations. 

Here k is the static dielectric constant and E is the energy of the donor 
state relative to the conduction band minimum. In silicon, for example, 
the function F corresponding to the minimum on the (0,0,1) axis satisfies 
the equation 


2mi dz^ 


2m, dy^/ kv^ 


F(t) = EF(t), 


(5.4) 


[See Eqs. (4.1) and (4.2).] 

It is convenient to adopt the following normalizations: 

I f k,(r)|*dr = 1, (5.5) 

Jcell 

j |F(r)|»dr = 1, (5.6) 

where is the volume per unit cell of the crystal. 

Before going on to the derivation of the wave function (5.1), let us 
give a qualitative picture of what it means. 

Equation (5.4) evidently is very similar to the wave equation for the 
hydrogen atom, the difference being that the free electron mass mo is 
replaced by two different masses, m, and m,, and has been replaced by 
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the much smaller quantity e^/x. As a result the lowest solution of Eq. 
(5.4) has a mean ‘‘Bohr radius of the order of 

^ (wVw^o) (^-7) 

where m* is an appropriate average between mj and mj and ao is the nor- 
mal Bohr radius. Since m*/mo is substantially less than 1 in both mate- 
rials [see Eqs. (4.2) and (4.3)] and k has the values®® 


Ksi = 12.0 (5.8) 

ICGe = 16.0 (5.9) 

reasonable choices of m* lead to the following magnitudes for a*: 

Osi* « 20 A (5.10) 

OGe* «45A. (5.11) 


These are enormous orbits, extending over several thousand crystal cells. 
Because of the heavier longitudinal mass, the wave functions Fj(r) 
do not have spherical symmetry but are somewhat compressed in the 
direction associated with the heavy mass. This compressive effect is 
however much smaller than the mass ratio mi/rrit might suggest [see 
Eqs. (6.5)]. 

The Bloch waves ^;(r) have the well-known form 

(5.12) 

in which fcy is the position of the jth minimum and Uj{r) is a completely 
periodic function. The functions Uj{t) vary quite rapidly within each cell, 
particularly near the nuclei where they resemble atomic wave functions. 
If we now consider one term 

x,(r) = f,(r)^,(r) (5.13) 

in the impurity wave function (5.1), we see that it can be well approxi- 
mated by 

Fj(Ta)<Pj(T) (5.14) 

in a region consisting of just a few lattice cells. Here ro is some point within 
this region. This shows that locally an electron described by (5.13) moves 
as though it were in the Bloch state ^i(r). However, if the electron is 
followed over distances of the order of a*, its wave function is slowly 
modulated by the function Fj(r) as a result of the weak Coulomb attrac- 
tion of the donor ion. In brief, (5.13) can be visualized as the short wave- 


H. B. Briggs, Phys. Rev, 77, 287 (1950). 
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length pattern, ^;(r), modulated in amplitude by the long wavelength 
envelope function Fj{x). (See Fig. 10.) 

Thus far we have discussed only one term of the sum (5.1). We shall 
show shortly that in fact any one term of this type represents an approxi- 
mate solution of the Schroedinger equation with some energy E. As a re- 
sult of the equivalence of the N conduction band minima, there obviously 
will be similar solutions, having the same energy E, associated with the 
other minima. The actual solutions of the impurity state problem are 



Fig. 10. Schematic representation of the donor wave function Xj(r)* 

given approximately by linear combinations of the degenerate functions 
which have the correct symmetry properties. For example, the ground- 
state wave function has the complete tetrahedral symmetry of the Hamil- 
tonian of the impurity state. This means that the coefficients ay in (5.1) 
are all equal if the phases of and F,(r) are chosen in a natural way.®^ 
We now turn to the derivation of Eq. (5,1) for the impurity state 
wave function. 

a. Single Band Minimum at\s. = 0 

Equations of the type (5.3) were first obtained by Wannier®^ with the 
aid of the functions which have become known generally as Wannier 
functions. A treatment of the impurity state problem along these lines 
may be found in a recent article by Slater,®® Here we shall base the discus- 
sion on the more familiar Bloch waves. 

For simplicity we consider first a semiconductor whose conduction 

Fy(0) and uy(0) real and positive for all j. 

” G. Wannier, Phys. Reo. 52, 191 (1937). 

” J. C. Slater, in ‘‘Encyclopedia of Physics*’ (S. Flugge, ed.), Vol, XIX, p. 1. Springer, 
Berlin, 1956. See especially Sections 25 and 26. 
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band has a single minimum at k = 0. We shall use the reduced zone 
scheme in which k is confined to the first Brillouin zone and the various 
energy bands are denoted by an index n. The lowest conduction band 
will be labeled by n = 0, and we shall take its minimum as our zero of 
energy. 

In the absence of an}^ impurities, the Hamiltonian for a conduction 
electron is 

ill which F(r) is an effective periodic potential. The eigenfunctions are 
the Bloch waves 

lAn,k(r) = ^ M„,k(r)e*-' (5.16) 

which we shall take to be normalized to unity in a large cube of volume 
V = The Wn,k(r) are periodic functions. We denote the energy corre- 
sponding to by £n(k). 

If now one of the atoms of the perfect lattice is replaced by a positive 
singly charged impurity ion, the Hamiltonian for an extra electron will 
be 

H = H^+ U{r) (5.17) 

where, at large distances from the impurity ion, one can write 

(.’(r) = - (5.18) 

K being the static dielectric constant. The reason for using the static 
dielectric constant is that the frequency of the impurity electron will 
turn out to be much lower than that of other electrons in the crystal, 
whose polarization gives rise to k in a homopolar material. The true 
U (r) will be a rather complicated potential in the immediate vicinity of 
the impurity ion, but we shall suppose, and to some extent verify later, 
that this region does not play an important role and shall use (5.18) for 
all values of r. 

We consider now the motion of the impurity electron as governed by 
the Hamiltonian H in (5.17). We expand its wave function yp in terms of 
the eigenfunctions ^„,k of the unperturbed Hamiltonian Hq: 

^ ^ 4n(k)^„,k. (5.19) 

For a discussion of the role of the dielectric constant in this problem, from the point 
of view of the entire many electron system, see W. Kohn, Phys. Rev, 106 , 509 (1957). 
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Here k runs over all points inside the first Brillouin zone which are con- 
sistent with periodic boundary conditions imposed on the surface of our 
large cube: 

kjj = 2Trla 


la being an integer. Substituting (5.18) into the Schroedinger equation 


{Hq + f/)^ = Ef!/ 

and taking the scalar product with ^n,k> we obtain 

(i?„(k) - £)^„(k) + 2 (nk|f/in'kO^»'(k') = 0 

n'.k' 

where 

(nkit/|»'k') = j 

= y I (- dr. 


(5.20) 

(5.21) 

( 5 . 22 ) 


Let us study this matrix element a little more closely. We note first that 
w*n,kWn'.k' is a periodic function and can be written as a Fourier series 
with the use of the reciprocal lattice vectors K,: 

W*n.k«n'.k' = ^ (5.23) 

p 

which gives 

(„k|f/|«f) - ^ (- tt’ - ti- )' <*■“> 


From the orthonormality of the ^n.k we find by integrating (5.23) under 
the condition k = k' = 0, 


C®n,0;n',0 


1 n = n' 

0 n ^ 


(5.25) 


It is clear that if the attractive force of the impurity ion is imagined 
to become weaker and weaker, the wave function ^ will more and more re- 
semble a Bloch wave near the conduction band minimum. Since the actual 
Coulomb force is strongly reduced by the dielectric constant, we antici- 
pate that the system of Eqs. (5.21) will have solutions in which u4n(k) is 
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negligible unless n = 0 and k is small. This suggests that in Eqs. (5,21) 
we display the terms with n = 0 in a prominent way: 

(Eo(k) - E)^o(k) + y (Ok|[/|Ok')Ao(k') 

F 

+ ^ (Okl[/|n'k')A„<(k') = 0 (5.26a) 

k' 


(^„(k) - B)A»(k) + y (nklC;|Ok')Ao(k') 

V 

+ ^ (nk|t7|»'k')^»'(k0 = 0 (ra 0). (5.26b) 

nVO 

k' 


As a first approximation let us set An(k) = 0 for n ^ 0 and replace 
EqOs) and (OkjJ/lOk) by their expressions for small k and k': 

^o(k)^^F (5.27) 

(0k|f/|0k') (5.28) 

Substituting into (5.26a) we find 


This will be recognized as the Schroedinger equation in momentum 
space of an electron of mass m* moving in the attractive Coulomb poten- 
tial — b^/kv. We note, however, that in (5,29) k and k' are restricted to lie 
within the first Brillouin zone. If this restriction is ignored, the solutions 
of (5.29) are almost zero for k outside the first Brillouin zone in any case, 
so that the restriction is not of great importance. For the time being, 
therefore, we replace (5.29) by the equation 




, Ao(k') = 0. (5.30) 


To transform this equation into a Schroedinger equation in coordinate 
space we introduce the Fourier transform of Ao in the usual way, namely, 

Fir) = p y Ao(k)e*'. 

all k 


(5.31) 
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Multiplying (5.30) by e**' and summing over k, we obtain 

a hydrogen-like so-called effective mass equation. The energy spectrum 
of (5.32) is of course 

no^ 2h^m*’ “ 1, 2, ... . (5.33) 


The eigenfunctions are familiar from the theory of the hydrogen atom. 
For example, the normalized fimction for the ground state is 


F(r) = 


1 

(7ra*8)i 


g-r/a* 


where the “effective Bohr radius '' a* is 


a 


* 


h^K 

m*e^ 


(5.34) 


(5.35) 


The total wave function of the impurity electron is, of course, not 
just F(r), but rather, in the present approximation 

^(r) = ^ A„(k)^„,k(r) (5.36) 

k 

^ X ^o(k)Mo.k(r)e*-^ 

k 

[see Eq. (5.18)]. Anticipating again that Ao(k) will be confined to the 
vicinity of k = 0, we may write 

Wo»k “ Wo,0 “f“ ka(dWo,k/^^a) (5.37) 

in (5.36) and, to a lowest approximation, neglect all but the leading term 
Uo,o> The result is 

^(r) « ^ Mo.o(r) ^ ^„(k)e*-^ (5.38) 

k 

= Wo.o(r)F(r) 

which is the form we wanted to derive [see Eq. (5.13)]. 

We now want to verify the two presumptions on which our derivation 
was based, namely that Ao(k) extends over only a small fraction of the 
first Brillouin zone and that the neglected amplitudes An(k), n 0 are 
in fact negligible. 
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If we invert Eq. (5.31) and use the solution (5.34) for F(r), we find 


Ao{k) 


Stt* 1 

+ (l/a*)2]2 


(5.39) 


so that Ao(k) extends appreciably in momentum space only to values of 
k less than or equal to 

k = 1/a*. (5.40) 


This must be compared to the mean value of k at the zone boundary, 
obtained by replacing the first Brillouin zone by a sphere of the same 
volume 


kzB = ^24iTr^/a = 6.2/a 


(5.41) 


where a is the lattice parameter. In silicon and germanium, a* has the 
approximate values 20 and 45 A, respectively, whereas a has the values 
5.42 and 5.62 A. Therefore, 


k/kzB — 


0.04 in silicon 
0.02 in germanium. 


(5.42) 


At the zone boundary Ao(k) has dropped to the fraction {k/kzsY of its 
maximum value, which is quite negligible for both materials. 

An estimate of An(k) for n 0 can be made from Eq. (5,26b). Treat- 
ing An(k) as small quantity, we have 

A„(k) = - _ -E Y( (»»k|f:|0k')A„(k') (5.43) 


which gives, in order of magnitude, 


l^nl 


AEa* ' "‘I* 


(5.44) 


Here Eo and AE are the ionization energy and band gap energy respec- 
tively. Thus in the present theory the ratio iAn|/|Aoj is a very small 
quantity, of the order of 10^® to 10”^, in silicon and germanium. 


b. Several Conduction Band Minima 

We mentioned in Section 4 that in both silicon and germanium the 
conduction band has several minima located at equivalent points in 
k space. This requires some relatively minor modifications of the simple 
effective mass theory which has just been developed.®®"®® 

« C. Kittel and A. H. Mitchell, Phys, Rev. 96, 1488 (1954). 

M. Lampert, Phya. Rev. 97, 352 (1955). 

37 W. Kohn and J. M. Luttinger, Phys. Rev. 97, 1721 (1955); ibid. 98, 915 (1955). 

3« W. Kleiner, Phys. Rev. 97, 1722 (1955). 
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In the present case we can find approximate solutions of the Schroed- 
inger equation consisting of Bloch waves associated with the vicinity 
of any one of the conduction band minima, say the one located at ky. 
We find, in complete analogy with the derivation leading to Eq. (5.38), 
that such solutions are given approximately by 

Xj(r) = Fj{x)ipj{x) (5.45) 

where <pj is the Bloch wave at the conduction band minimum k, and 
Fj(r) is an envelope function satisfying the effective mass equation 

E>(k “ ky) is the energy of the conduction band near k = ky. 

When two effective masses are identical, as in silicon and germanium, 
Eq. (5.46) has axial symmetry and its solutions can be labeled by a 
magnetic quantum number m. Solutions corresponding to different values 
of I ml are no longer degenerate as in the case of an isotropic mass, but the 
twofold degeneracy arising from + and — m remains. 

Consider now, for example, the lowest energy solution of (5.46), which 
has m = 0. Clearly there also are (iV — 1) equivalent solutions, asso- 
ciated with the {N — 1) other conduction band minima, all of which 
have the same energy in the effective mass approximation. This N-fo\d 
degeneracy is reduced, however, when corrections to the effective mass 
formalism are taken into account. The degeneracy which remains can 
be ascertained solely from the symmetry of the true Hamiltonian of the 
impurity state — invariance under the tetrahedral group — without 
detailed knowledge of the nature of the corrections. This may be achieved 
by the following well-known method. The N degenerate functions clearly 
form the basis for a representation R of the group Tj. This representation 
may be “reduced'^ into the irreducible representations Ai, A 2, Ej Ti, 
and T 2 of the group Td,^^ 

R = UiAi “b n^A^ H” nzE -f- n^T\ -|- /15T2 (5.47) 

where the n, are integers. Let us call d* the dimensionality of the represen- 
tation associated with ni. Regardless of the detailed nature of the cor- 
rections to the effective mass formalism, the N originally degenerate 
levels will then split into ui levels with degeneracy di, with degeneracy 
d2, etc. 

Let us illustrate these remarks by using silicon, for which W = 6, 

H. Eyring, J. E. Walter, and G. E. Kimball, “Quantum Chemistry,” p. 388. Wiley, 
New York, 1944. 
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as an example. For definiteness, let us take the arbitrary phase factors 
of the xy in such a way that all Xj have the same value at r = 0.®^ We 
write the linear combinations corresponding to the various irreducible 
representations in the form 

6 

^ z = 1, . . . 6 (5.48) 


where j = 1, . . . 6 are taken to correspond to the minima (fco,0,0), 
( — fco,0,0), . , . (0,0, — A:o), respectively. Elementary group theory then 
yields the following linear combinations corresponding to the different 
irreducible representations : 


a/>’ 

1/V6(1, 1,1, 1.1,1) 

(^i) 


Kl, 1,-1, -1,0,0) 1 

{E) 

a/’) 

^(1,1,0,0,-1,-1) J 


l/\/2(l,-l,0,0,0,0) ] 



1/V'2(0,0,1,- 1,0,0) 

■{T,) 


1/V2(0,0,0,0,1,-1). J 



(5.49) 


No tetrahedrally symmetric perturbation can lift the double degeneracy 
of and or the triple degeneracy of and Therefore, 

we shall in reality have a singly, a doubly, and a triply degenerate level 
associated with the lowest eigenvalue of the effective mass equation. 
The spacing between the three levels will, of course, be the smaller, the 
better the effective mass theory accounts for the true state of affairs. 

The situation is similar for the higher eigenvalues of the effective 
mass equation. Further details can be found in the article of reference 37. 

c. Matrix Elements for Radiative Transitions 

We have seen in Section lb, that a great deal of information has been 
obtained from infrared absorption experiments in which a donor electron 
makes a transition from the ground state to an excited state. The transi- 
tion rates are proportional to the squares of the matrix elements of mo- 
mentum which we shall now evaluate. 

We begin with the case of a simple band having a minimum at k = 0, 
in the vicinity of which the energy is given by 

Eik) (5.50) 

Let us denote the ground state by and the excited state of interest 
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by To calculate the momentum matrix elements it is essential to 
use the expressions 

^ A«("ki(»o.k(r), I = 0 or e (5.5V) 

k 

for the wave functions rather than the less accurate expressions of the 
form F(^^(r)wo.o(r).^® Making use of the fact that Pa = {l/i)d/dxa is a 
periodic operator we find, 

k V 

= ^ Ao<''>*(k)Ao''Hk)(i;'o.k,p«iAo,k) 

k 

= ^Yi 

k 

= 5 (/’«>(r),p«F«(r)). (5.52) 

In view of the commutation relationship 

[xa,H\ = = — Pa (5.53) 

in which H is the total Hamiltonian appearing in (5.20), we can obtain 
the alternative result 


(^(•),p.^(.)) = ^ (^(«),[Xa,H]^‘*>) 

= ^ (£<•> - fi<«)(^(w,a:«,^<*)) 

vn 

= g (5.54) 

Here and are the energies of and respectively. 

It should be noted that neither (5.53) nor (5.54) require a knowledge 
of the Bloch wave at k = 0. 

In the case in which there are N equivalent conduction band minima 


If the latter type of wave functions are used one obtains instead of (5.52) the 
result paF^*^(t))] i.e., one loses the essential factor m/m*. 
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we may let the ground and excited states have the wave functions 
and of the form 

N 

^(0 = 2 (5.55) 

y-i 

A simple extension of the derivation just given leads to the result 

N 

g (£:C5 - £<»)) ^ (5.56) 

t=\ 

6. Numerical Results and Comparison with Experiments 

We shall now compare the predictions of the simple effective mass 
theory, developed in the preceding section, with experiment. 


a. Energy Levels 

To obtain the energy levels in the effective mass approximation we 
must determine the eigenvalues of Eq. (5.4), 


2mi dz^ 


2mt 



F{r) = EF(t) 


(5.4) 


in which the parameters mi, mt and k have the values given in Eqs. (4.2), 
(4.3), (5.8), and (5.9). These values have been determined by entirely 
independent experiments and will not be regarded as “adjustable.'’ 

Since (5.4) has axial symmetry it can be partly separated and reduced 
to a two-dimensional partial differential equation. The latter must then 
be solved approximately, for example by variational methods.*®”’® The 
solutions are, of course, rather different from simple hydrogenic ones 
because of the large ratio of mi/m* in both silicon and germanium. Never- 
theless, it is convenient to label them by the familiar atomic notation 
(Is, etc.) which indicate the hydrogenic functions into which they go in 
the limit where > 1. The magnetic quantum number m refers, 

of course, to the remaining axis of symmetry. The results are given in 
Table V. Let us compare them with the experimental data. 

In germanium only the ionization energy of the Is state has been 
measured. For the different impurities, it lies between 0.0096 and 0.0127 
ev (see Table 1). This is in rather good agreement with the effective mass 
value of 0.0092 ev. 

In silicon the effective mass value of 0.029 ev for the ground state is 
in good agreement with the ionization energy of Li (0.033 ev), but agrees 
poorly with those of P, As, Sb, and Bi (0.045 to 0.067 ev). Since the orbits 
in silicon are smaller than those in germanium, it is not surprising that 
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Table V. Donob Spectra in Effective Mass Approximation 


Energies in units of 0.01 ev® 


State Silicon 


Germanium 


Is 



-2.9 

± 

0.1 (37) 

-0.92 

± 

0.02 (42) 

2p, 

m 

= 0 

-1.09 

± 

0.02 (41) 

-0,45 

± 

0,02 (42) 

2s 



-0.88 

± 

0 06 (37) 





m 

= ±1 

-0.59 

± 

0.01 (41) 

-0.160 

± 

0.003 (42) 

3p, 

m 

= 0 

-0.57 

+ 

0.06 (41) 

-0.235 

± 

0.02 (42) 

3p, 

m 

= ±1 

-0.29 

± 

0.005 (41) 

-0.085 

± 

0.005(42) 


® Numbers in parentheses refer to footnote references in text. 

the effective mass theory is less successful. We shall return to this com- 
parative failure of the theory in the following Section 7. 

Turning now to the excited states in silicon (see Table II), we note 
the following experimental values of the differences in energy between 
successive excited states : 


P 

As 

Sb 


0.050 

0.053 

0.047 

ev 

0,031 

0,032 

0.034 

ev 

0.020 

0.023 


ev 


It is encouraging to observe that the first two differences in energy are 
very similar for all donors whose excited states have been studied. This 
would, of course, be expected from the effective mass theory. Before 
comparing the experimental results with the calculations in more detail, 
we first note that the matrix elements for radiative transitions from the 
ground state vanish in the effective mass approximation unless the en- 
velope function of the excited state in question has odd parity [see Eq. 
(5.56)]. Approximate calculations of the total intensity give the results 
listed in Table VI.^^ Only limited reliance can be placed on these results. 

Table VI. Calculated Optical Absorption Intensities in Siucon 


Transition 

Relative intensity 

Is — ► (2p, m = 0) 

4.0 

Is — > (2p, m = ±1) 

10.6 

Is — > (3p, w = 0) 

0.4 

Is— ► (3p, m == ±1) 

3.1 


W. Kohn, Phya, Rev. 98 , 1856 (1955). 
W. Kohn, unpublished. 
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for the matrix elements involve the Is function, which, as we have seen 
in connection with the Is energy, is understood only imperfectly. The 
small intensity Is (3p, m = 0) is probably significant however. More- 
over, since the (3p, m = 0) level should be almost coincident with the 
‘‘strong'^ (2p, m = ±1) level (see Table V), it is quite likely that it 
would escape detection. Consequently, we are led to compare the first 
two experimental spacings (6,1) with the following theoretical energy 
differences : 

(2p, m = ± 1) — (2p, m = 0) lO. 060 ev 

(3p, m = ±1) — (2p, m = ±1); 0.030 ev. 

The agreement is strikingly good and is one of the brightest points in the 
field of impurity states. The comparison between experiment and theory 
is also shown in Fig. 3. 

It is easy to present two good reasons why the effective mass theory 
should apply better to the excited states than to the ground state. First, 
their orbits are larger; second, their envelope functions have odd parity 
and thus vanish at the donor nucleus. For both these reasons the orbits of 
the excited state penetrate quite negligibly into the region near the donor 
ion where our expression — e^//cr for the potential energy presumably 
breaks down. 


6. Wave Functions 

The wave function of greatest interest is that of the ground state. 
It has the form 

N 

^(r) = J «,f,(r)^,(r) (6.3) 

[see Eq. (5.1)]. 

Let us begin by discussing the coefficients a^. As we have seen, the 
ground state is iV'-fold degenerate in the effective mass theory, correspond- 
ing to N different sets of ay. This degeneracy is partly lifted when cor- 
rections to the effective mass formalism are taken into account. However, 
it is not certain on a priori grounds to which representation the resulting 
lowest level belongs. One would be inclined to expect that the completely 
symmetrical, nondegenerate state Ai [see Eq. (5.49)], in which all ay are 
equal has the lowest energy. Fortunately there is a crucial test of this 
assumption. We can see from Eq. (5.49) that only the state Ai has a 
nonvanishing wave function at the donor nucleus, which produces a 
large hyperfine interaction with the spin of the donor electron. In silicon, 
the experimental observation of this interaction, which has been men- 
tioned in Section 2a and will be discussed further below, provides con- 
clusive proof that the ground state has the complete tetrahedral sym- 
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metry characteristic of A\. Donor spin resonance has not been observed 
in germanium, although it has been looked for^® vigorously. Hence, no 
clear-cut decision can be made in this case. One might wonder if the 
assumption of a degenerate ground state belonging to a multi-dimensional 
representation of Td would account for the failure of spin resonance ex- 
periments in germanium. 

We now turn to the envelope functions Fy We consider any one of the 
conduction band minima and take the direction defined by its k vector 
as the z axis. To a very good approximation the corresponding normalized 
wave function can be represented by 


p = —4— (6.4) 

(Tra^o)* 

This functional form, which is of course exact in the limit of equal effec- 
tive masses, was also found to be very accurate in the limit mi! rat —> 00 .*^ 
A variational calculation gives the following values for a and b : 

Osi = 25.0 X 10“® cm, b&i = 14.2 X 10“® cm , . 

Og* = 64.5 X 10“« cm, bo. = 22.7 X 10-» cm. 


We see that the wave functions extend further in the x and y directions, 
which correspond to the smaller transverse mass, than in the z direction, 
which corresponds to the larger longitudinal mass; however, the distortion 
is relatively small. Even in germanium, which has a mass ratio of 20 to 1, 
the ratio of a to 6 is only 3 to 1. Naturally we obtain a tetrahedrally 
symmetrical structure when the full wave function (6.3) is considered, 
which may be pictured as approximately spherical with a “radius” given 
roughly by the larger length o. 

Finally we come to the Bloch wave functions <py They have already 
been discussed in Section 4. We want to add only one important item of 
information. By studying the rate at which partially aligned Si^® nuclei 
relax by means of the hyperfine interaction with conduction electrons, 
Shulman and Wyluda^^ obtained an estimate of the square of the wave 
functions associated with silicon conduction band at the positions of the 
Si nuclei. Their results are best expressed in terms of a dimensionless 
quantity which measures the extent to which the conduction band wave 
functions are concentrated at the lattice sites r* : 


V — 


ky(r)l^Av 


= 186 ± 18. 


( 6 . 6 ) 


We are now ready to compare the theoretical value of at the posi- 


G. Feher, private communication. 

** R. G. Shulman and B. J. Wyluda, Phys. Rev. 103, 1127 (1956). 
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tion of a P donor in silicon, with the experimental value obtained from 
the hyperfine splitting of the donor electron resonance (Section 2a). 
Using (6.3) with all = \/^7 closing a, ft, and i? as given by Eqs. 
(6.5) and (6.6), we obtain 

\m\^ = 6lFX0)l^|,.y(0)l2 

= 0.042 X 10^^ cm-®. (6.7) 

This is much smaller than the experimental value, namely 0.44 X 10^^ 
cm“* (Table III). The difference represents a serious failure of the effec- 
tive mass theory. We shall see in the following section how most of this 
discrepancy can be removed by a simple correction to the formalism. 

7. Corrections to the Effective Mass Formalism 

a. General Consider atims 

We have seen some of the successes and failures of the effective mass 
formalism in the preceding section. They are roughly what we might have 
expected. For states having large orbits, which do not penetrate into the 
immediate vicinity of the donor ion appreciably, the theory agrees well 
with experiment; for states having smaller orbits and substantial ampli- 
tudes near the donor ion, notably the ground states in silicon, the agree- 
ment is poor. 

What are the chief reasons for the failures of the effective mass for- 
malism? We have based our discussions on the Hamiltonian 

[see Eqs, (5.17), (5.15), (5.18)] and have shown that it leads to an effec- 
tive Hamiltonian, 

(7.2) 

The reason for concentrating on P (Z = 15) is that there exists a mathematical 
limit where the effective mass wave function (6.3) should become exact everywhere, 
even at the position of the P nucleus. This is the limit in which the charge on the 
P nucleus is imagined to be only inffnitesimally larger than that on the Si nuclei 
(i.e., Zp « 14-1-5). Ill limit the orbits become very large and the wave function 
is, locally, just like a superposition of silicon conduction band functions. On the 
other hand, consider another donor such as As (Z = 33). Here even in the limit in 
which the Coulomb attraction is imagined to be very weak (Z = 32 -h 5), the donor 
wave function in the cell of the As ion will resemble conduction band wave func- 
tions of Ge (Z = 32) more than those of Si. This fact is not allowed for in the effec- 
tive mass expression (6.3). 
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[see Eq. (5.46)]. The errors involved in passing from (7.1) to (7.2) were 
discussed in Section 5, and are very small indeed for silicon and ger- 
manium, being almost certainly less than 1 %. The trouble must therefore 
lie with the original Hamiltonian (7.1), Unfortunately a thorough analysis 
of these difficulties has not been published yet; therefore, we must be 
content with some rather incomplete remarks. 

The first question which may be asked is whether the impurity elec- 
tron can be described adequately by any single particle Hamiltonian. 
After all, it interacts constantly and strongly with a large number of 
electrons which lie within its orbit. It is by no means obvious that the 
actual many-particle problem can be replaced adequately by a single 
particle equation. A partial answer to this question was given by the 
author,^® who showed that the entire many-electron Hamiltonian can be 
reduced to a single electron effective mass Hamiltonian of the type H^tt 
of Eq. (7.2) in the hypothetical case in which the charge on the impurity 
nucleus exceeds that of the other crystal nuclei by an infinitesimal 
amount. The same result can also be derived^’ for sufficiently highly 
excited states in the case in which the donor nucleus possesses a finite 
excess charge. This is a satisfactory result in view of the good agreement 
between the effective mass theory and experiments involving the excited 
donor states. 

For the ground state of the donor, the situation near the donor ion 
is quite complicated. However, for r '> r, [(47r/3)r/ = atomic volume] it 
can be shown^^ that F(r) satisfies the effective mass equation 

^ “ S] ^ 

provided only that the fractional change of F{t) over a lattice spacing 
is small. Here E is the actual (observed) energy. It is to be expected that 
(7.3) will break down for r ^ r,. For example, the concept of the dielectric 
constant is meaningful only for sufficiently large distances. Moreover, 
the donor ion may distort the lattice significantly in its immediate 
vicinity. 

6. Corrected Wave Functions^^ 

For definiteness let us now consider a P donor in silicon. On the basis 
of Eq. (7.2) we can make the following simple correction to the effective 

W. Kohn ref. 34. 

Unpublished calculations by the author. 

** The following discussion is based on W. Kohn and J. M. Luttinger, Phya. Rev. 97, 

883 (1955). 
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mass wave function in the exterior region r r,: Instead of solving (7.3) 
over all space with the boundary conditions F{<x>) =0 and F(0) = finite, 
which determine the eigenvalue we solve (7.3) in the exterior region 
using only the one boundary condition F( oo ) =0 and the observed energy 
Eoh»* For simplicity let us also replace rrti and by a single mean mass 
m*. The latter may for example, be chosen so as to give the same energy, 
in the effective mass approximation, as mi and We then have in 
the exterior region 


6 


r^Tr.: = ~F(t) vAr) 

(7.4) 

j=l 

(7.5) 

f(oo) = 0. 

(7.6) 


Now £ob» (—0.045 ev) is substantially lower than the effective mass 
energy (—0.029 ev). Therefore, the solution of (7.5) and (7,6) falls off 
more rapidly than the effective mass wave function for large r, and since 
J5obs is not an eigenvalue of the effective mass equation, F(r) approaches 
infinity as r becomes small. Both these effects tend to concentrate the 
donor electron near the donor ion much more than in the effective mass 
formalism. 

The asymptotic behavior of F(t) for large values of r is 

Fir) ^ (7.7) 

where b* is related to the ‘‘Bohr radius’^ a* associated with the effective 
mtiss approximation by the equation 

6* = a*(E,um^/Eo^,)K (7.8) 

We shall treat the interior region, r ^ r*, rather crudely. Let us write 
the Schroedinger equation for ^ in the form 

(- ^ + F + [/) ^ = (7.9) 

where V is the effective periodic potential for the Bloch electrons near 
the conduction band minima and U is the perturbation caused by the 
extra charge on the P nucleus. Energies are measured from the bottom 

Other ways of choosing m*, for example by the condition that ^"(0) should have the 
same value as that calculated with mi and mt give very nearly the same result. 
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of the conduction band. If this extra charge is imagined to become 
infinitesimal, (7.9) clearly reduces to 

(-Si 

and the solution of (7.10) with the correct tetrahedral symmetry is 

6 

^ = C J ^,(r). (7.11) 

J = 1 

With the actual extra charge e on the P nucleus, U and £?ob. are still small 
compared to V, Thus we shall write the solution of (7.9) in the form 

6 

^ = F(r) ^ vm (7.12) 

where F(r) is not expected to vary appreciably over the distance r,. 
Rough estimates*®’^* indicate a variation of probably no more than 
10-20%. Although (7.12) still has the form of an effective mass wave 
function, it should be emphasized that we should not expect F{r) to 
satisfy an effective mass equation for small r. As we lack detailed knowl- 
edge of F(r) for r ^ r„ we can only set 

r<r.: F(r) « F(r.) (7.13) 

where |r,| = r*. 

It may appear that there is a contradiction between this picture of a 
slowly” varying F(r) and the statement following Eq. (7.6) that F{t) 
approaches infinity for small r. Actually there is no difficulty. The in- 
crease of F{t) as r — > r, from the outside is quite rapid on the scale of a*, 
the ^‘Bohr radius” of the effective mass equation, but is slow over dimen- 
sions of the order of r,. 

A schematic picture of the envelope function F(r), corrected for devia- 
tions from the effective mass formalism, is shown in Fig. 11. 

For other donors in silicon, such as Li, As, Sb, and Bi, the discussion 
of the exterior region is the same as for P. One merely has to use the ob- 
served energy of the ground state for each donor. This makes F(r,) a 
monotonic function of Eob.- The passage into the interior region is much 
more difficult, however, for this region is now occupied by an ion which is 
completely different from the atoms of the host crystal and the form of 
the donor wave function is quite different from the expression (7.12). 
An extremely qualitative treatment of the interior region for these donors 
is given in ref. 48. 



SHALLOW IMPURITY STATES IN Si AND Ge 


293 



Fig. 11. Comparison of the corrected envelope function and the effective mass 
envelope function for a P donor (schematic). 


CONDUCTION 
BAND EDGE. 
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-4*5 (experimental)" 


Fig. 12. Theoretical spectrum of a P donor. The numbers in parentheses refer to 
number of (approximately) degenerate states, spin degeneracy not included. 
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We remarked in Section 5b, that in the effective mass approximation 
there are 6 degenerate wave functions corresponding to the lowest energy 
eigenvalue. Now we have seen that the completely symmetrical state 
(belonging to representation A i) has in fact an energy that differs quite 
appreciably from the effective mass value. The other 5 states have wave 
functions which vanish at the donor nucleus. One would, therefore, ex- 
pect them to be less affected by the breakdown of the effective mass 
formalism near this nucleus. Their energies should then be more nearly 
equal to the effective mass energy. 

Similarly, of the six ^‘2s^^-like states only one is expected to be de- 
pressed appreciably, by an amount which can be calculated from the 
observed depression by means of the quantum defect method. (See 
ref. 48.) The resulting spectrum for a P donor is shown in Fig. 12. 


c. Comparison with Experiment 

In Section 6b we calculated |^(0)|2 for a P-donor electron by means 
of the effective mass formalism and found it to be 10 times smaller than 
the experimental value. We have also noted in Section 6a that the ob- 
served ionization energy is about 50% larger than that calculated by the 
effective mass theory. We shall now see that the first large discrepancy 
disappears when the experimental ionization energy, rather than the 
theoretical value, is properly included. 

Following the discussion of the preceding subsection we have 


1^(0) = 


6|F(0)|2l^y(0)|2 

6iF(r.)121^y(0)l^ 


(7.14) 


Now when the properly normalized®^ F(r) is calculated from Eqs. (7,5), 
(7.6) we find that F(r*)^ is 9.8 times as large as the value of |F(0)|^ cal- 
culated in the effective mass approximation. This large enhancement is 
produced by the relatively small shift in energy because F(r) tends to 
become infinite for small r. When this correction is applied to the effective 
mass result (6.7) one obtains 


|^(0)|2corrected = 9.8|^(0) | = 0,41 X 10^^ cm“®. (7.15) 

This agrees very well with the experimental result, namely 0.44 X 10^^ 


^0 For normalization purposes one can use the singular solution of (7.5), (7.6) for all r, 
since it is square integrable and the interior’' region makes a negligible contribu- 
tion because of its small volume. 
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cin“®, in fact somewhat better than one has reason to expect. For example, 
the choice of r, as the dividing line between interior and exterior regions 
is somewhat arbitrary. Nevertheless, it is gratifying that this major 
failure of the effective mass equation can be remedied so simply and 
naturally. 

For other donors the very rough discussion of 1^(0) given in ref. 48 
is found to be in reasonable agreement with experiment. Let us take the 
case of Li for example. Here the ionization energy is very close to the 
effective mass value (0.033 and 0.029 ev, respectively) so that the “en- 
hancement factor’^ 9.8 of Eq. (7.15) is absent. Furthermore, Li is a very 
light atom so that the density of the donor electron is not as concentrated 
near its nucleus as near a silicon or phosphorus nucleus. Taking both 
these factors into account, the theoretical estimate for 1^(0) was 0.002 X 
10^^ cm""® which may be compared with the recently measured value of 
0.0032 X 1024 cm-2. 

Let us next ask how well the theory can account for the effects of 
the hyperfine interaction of the donor electron with the Si^® nuclei (see 
Section 2). For this purpose we can use the exterior form of the wave 
function 

6 

Hr) = ^ ^,(r) 

6 

= (7.16) 

3^1 

for all donors, where F(r) is of course computed with the ionization energy 
appropriate to the donor in question. The resonance frequency of a Si^® 
nucleus located at n is determined by j^(n)l® according to Eq. (2.2'). 
With the phase convention of footnote 39a, all six wy(r) are equal at each 
lattice site, say u{ti), so that we can remove the common factor w(n) 
from the sum in (7.16) and write 

6 

(7.17) 
y=i 

Here we see an interesting interference effect between the six wavepackets 
of which ^(r) is composed; The magnitude of yp{ti) depends not only on 
F{xi) and \u{ti)\ = [see Eq. (6.6)] but also on the extent to which the 
six terms in (7.17) interfere with each other. This interference is deter- 
mined by the magnitude fco of the six ky, that is, by the exact location of 
the minima of the conduction band along the (1,0,0) and equivalent axes. 


Hrd = -^F(ri)u(ri) ( V 


V6 
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Since fco is not yet well known, D. Schechter, in unpublished work, 
has plotted 


6 




(7.18) 


as function of ko, for lattice points located near the donor ion. The com- 
parison with the experimental results of Feher^® is complicated by the 
fact that one does not know, of course, which experimental resonance line 
corresponds to a specific fj. The general magnitudes of the experimental 
and theoretical values of are in good agreement, but a more de- 

tailed analysis is still in progress. Such an analysis might test not only 
the theoretical picture of the donor states but could also pin down the 
exact position fco of the minima of the conduction band. 

Let us turn next to the width of the electron spin resonance arising 
from the hyperfine interaction of the electron spin with the magnetic 
moments of the Si^® nuclei (see Section 2a). The shift of the resonance 
field of a particular donor electron arising from its interaction with the 
moments of Si^® within its orbit is given by 

A'H = (7.19) 


Here fi&i and Isi are the magnetic moment (—0.555 nuclear magneton) 
and spin (i) of Si*®, i runs over all lattice sites occupied by Si®® nuclei and 
nit is the magnetic quantum number of the Si*® nucleus at r^; A'H differs 
from donor to donor because of differences in r« and If one uses the 
fact that Tt and are random, one finds the expression 

= P (y Msi) (7.20) 

for the root-mean-square value of A'/f. Here / is the natural fractional 
abundance of the Si*® isotope (0.0468) and the sum now runs over all 
lattice sites except that of the donor. 

Equation (7.20), which again depends on fco, must be compared to 
the observed root-mean-square widths of the electron resonance lines. 
A preliminary comparison of this type has been reported®^ and appears 
promising. It was mentioned in Section 2a that the hyperfine splitting of 
the electron resonance arising from interaction with the donor nucleus 
is much larger (42 oersted for P) than the width of each multiplet («3 
oersted) associated with its interaction with ail Si*® nuclei in its orbit. 


W. Kohn, ref. 26. 
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This indicates a much weaker interaction with the Si^® nuclei. The differ- 
ence is due to the following two causes: first, the envelope function F(r) 
becomes exceptionally large at r = 0, as was mentioned in the discussion 
following Eq. (7.6); second while the planewave terms exp (tkyr^) are 
all in phase at u — 0, they will interfere substantially at most other 
lattice sites. It appears at present that the line widths of Sb, P, and A 
(Table IV) can be fitted within about 25% by a choice of ko in the vicinity 
of 0.7fcmax, Here km^. is the value of k at the zone boundary. Experimental 
and theoretical work on this question is still in progress and the quoted 
estimate of fco must be regarded as very tentative. 

IV, Structure of Acceptor States 

In this chapter we shall review our current understanding of the 
shallow acceptor states in silicon and germanium. While these states are 
analogous to the donor states in a general way, we do not have nearly 
as detailed a picture of them. The reason for this is twofold. In the first 
place the theory of the states involves the structure of the valence band 
near its top in an essential way and we shall see in Section 8 that this band 
is much more complicated and less well known than the structure of the 
conduction bands. Secondly, spin resonance experiments, which have shed 
so much light on the nature of donor states, have not been successful in 
the case of acceptor states. Fortunately the infrared absorption spectra of 
these states, which were first observed by Burstein and co-workers,^ have 
been determined recently by Hrostowski and Kaiser^* with remarkable 
precision. Their measurements promise to be an important aid in the 
further elucidation of the nature of the acceptor states. 

8. Valence Bands of Silicon and Germanium 

Acceptor states are wave packets made up largely of Bloch waves 
chosen from near the top of the valence band. To understand them one 
requires some information about the structure of this band. This informa- 
tion, which has for the most part been extracted from cyclotron resonance 
experiments, will now be summarized. A more detailed analysis of 
the valence bands of silicon and germanium can be found in articles by 
F. Herman, Dresselhaus et and E. Kane.®^ 

a, Silicon 

The valence band of silicon is shown schematically in Fig. 8. Its 
highest point is located at k = 0. If it were not for spin orbit coupling, 

G. Dresselhaus, A. F. Kip, and C. Kittel, Phys. Rev. 98 , 368 (1955). 

B. Lax, H. J. Zeiger, and R. N. Dexter, Physica 20, 818 (1954). 

K. O. Kane, /. Phys. Ckem. Solids 1, 83 (1956). 
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this point would be sixfold degenerate, including the double degeneracy 
arising from spin. The simplest way to understand this degeneracy is to 
consider the tight binding limit, in which the wave functions correspond- 
ing to the highest point go over into atomic 3p functions. The spin orbit 
coupling lifts the degeneracy partially and leads to the situation shown 
in Fig. 8. The top of the valence band remains at k = 0 and is fourfold 
degenerate, corresponding to atomic J = f states. Slightly below it is a 
twofold degenerate point corresponding to / = The energy splitting 
X between the two points has not yet been measured directly. Theoretical 
estimates give 

X « 0.035 ev (8,1) 


with an uncertainty of perhaps 20%. For small values of k, the energy E 
of the valence bands can no longer be expressed as a power series in kzj 
kyy and kg. Instead, it has been shown that it is a solution of a 6 X 6 
secular equation of the following form: 


Det \Di/^kak 0 - - ESjA = 0, j,j’ = 1, ... 6. (8.2) 


Here summation of a and jS over 1 to 3 is implied; is defined by 


I 0 i = 1, . . . 4 

I i = 5, 


(8.3) 


the coefficients 0 ^“^ are certain numerical constants characteristic of the 
material and of the dimension h^/vto; E is the energy reckoned from the 
top of the valence band. The 6X6 matrix D = Djj'^^kakfi is Hermitian 
and has the following structure;®* 



P/2 

R* 

S* 

0 

-iR*/V2 

iR*lV2 


R 

P/6 + 2Q/3 
0 

S* 

-(P-2Q)/3V2 
-i VI ft* 

iRly/2 

-i{P - 2Q)/3 

-i VI R* 

(P + Q)/3 


S 0 

0 s 

p/6 + 2Q/3 -ft 

-ft* P/2 

i VI ft _ iV2S 

i{P - 2Q)/3 V2 iR/y/2 

-iV2S 

V2 f V|ft 

-i(P - 2Q)3 V2 
-iR*lV2 
0 


0 (P + Q)/3 


(8.4) 


W, Shockley, Phys. Rev. 78, 173 (1950). 

J. M, Luttinger and W. Kohn, Phys. Rev. 97, 869 (1955), 
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where 

P(k) - - + (L + + 2MK\ 

7R 

Q(k) = ^ + V) + (8.5) 

ft(k) = -(*Ar/V3)(fcx — tky)k^j 

S{k) = (l/VT2)[(i- - M){k,^ - - 2t-iVfc,y. 

For silicon, the constants have the following values 
(L + 2M)/3 = (-5.04 ± 0.12) ^ 

(L - M)/3 = ±(L24 ± 0.46) ~ (8.6) 

N - (-7.49 ± 0.40)^. 

Unfortunately, the sign of (L — M)/3 is not determined by classical 
cyclotron resonance experiments. 

For k = 0, (8.2) has a fourfold root at = 0 and a twofold root at 
E = —X, as it should. For sufficiently small values of fc,®* it can be fac- 
tored approximately into the following 4X4 and 2X2 equations: 

Det \Djf^^kakfi - Edjf\ =0 j,/ = 1, . . . 4 (8.7a) 

and 

Det IDjr^^kakp - (E + \)Bjy\ = 0 i, f = 5, 6. (8.7b) 

Equation (8.7a) describes the vicinity of the very top of the valence band, 
the point B of Fig. 8, whereas Eq. (8.7b) describes the vicinity of the 
split-off point B\ 

The quartic equation (8.7a) can be factored and solved, giving 

E{k) = Ak^^[B%^ + + h^k/)]^ (8.8) 

each eigenvalue being doubly degenerate. The constants A, Bj C are 
related to L, Af, and N of Eq. (8.6) as follows: 

A-i(L + 2M) + ^ 

B = UL- M) (8.9) 

C2 = i[i\r2 ~ (L - My\ 

We see that E ^ k!^ in any given direction in k space, but since the curva- 

These values were kindly communicated by Dr. H. J. Zeiger of the M.I.T. Lincoln 
laboratories. 

More precisely, this means that E(k) — ^(0) must be small compared to X. 
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ture depends on the direction of k, the surfaces of constant energy are 
warped. 

b. Germanium 

The valence band structure of germanium is completely analogous to 
that of silicon. However the parameters are much more accurately known 
in germanium. They are:^^ 

L = (-32.0 ± 0,2)^ 

^ 2m 

M = (-5.3 + 0.2) (8.10) 

N - (-32.2 + 0.2) 


The sign ambiguity which existed in the case of silicon [see remark 
following Eq. (8.6)] has been resolved beyond doubt by band theoretical 
calculations. 

The spin orbit splitting at k = 0, X, is known from infrared absorption 
measurements to be 

X = 0.29 ± 0.02ev.®» (8.11) 


9. Effective Mass Equations for Acceptor States 


Because the valence bands of silicon and germanium are degenerate, 
the effective mass theory of acceptor states is appreciably more com- 
plicated than for donor states. Instead of a single partial differential 
equation, Eq. (5.4), one now finds a set of coupled partial differential 
equations. The reader is referred to papers by Kittel and MitchelP® and 
Luttinger and Kohn^® for derivations of these equations. They are 


" \i dX^) \i dX^ 


i) Fdr) - + ^F,(r) = EFj{i), 

J = 1, ■ • 


6. (9.1) 


Here the parameter cy, and X are the same as those occurring in the 

secular equation (8.2) for the determination of the energy of the valence 
bands, b^/kv represents the Coulomb interaction of the negatively charged 
acceptor ion and the hole, and E is the energy measured from the top of 
the valence band (E > 0 for bound states). 

The functions Ey(r) are envelope functions in the following sense. 
At k = 0 there are four degenerate Bloch waves . . . <pA, at the top 


“ A. H. Kahn, Phys. Rev. 97, 1647 (1955); W. C. Dash and R. Newman, Phys. Rev. 
99, 1151 (1955). 



SHALLOW IMPURITY STATES IN Si AND Ge 


301 


of the valence band, and two Bloch waves and at an energy X below 
the top, the latter being split off by spin orbit coupling. We take them to 
be normalized according to Eq, (5,5). The total wave function of an 
acceptor state is then given by 


Its normalization is 


6 

4>{t) = ^ 

j = l 


6 


/ l^(r)|^dr = 2 / \F,iT)\^dr. 

J-1 


(9.2) 

(9.3) 


The functions Fj{r) have the same general character as in the donor 
states. They extend smoothly over dimensions of the order of 15 A and 
more and can be thought of as generalizations of hydrogenic wave func- 
tions. On the other hand, the <Pj{r) are rapidly varying and completely peri- 
odic functions which have p character in the vicinity of any lattice site. 


10. Approximate Solutions and Comparison with Experiment 

Obtaining the Eq. (9,1) for the determination of the F{t) is a fairly 
simple matter. On the other hand, it is somewhat more difficult to find 
reasonably good solutions of these equations. Some exact solutions have 
been found, but only for special values of the parameters L, M, and N 
of Eqs. (8.5). In general, one has to rely on the variational method coupled 
with group theoretical considerations. A detailed analysis of the equations 
and their solutions will be given in a forthcoming thesis by Schechter.®® 
We shall merely sketch some of the salient features here. 


a. Germanium 


Let us begin with a discussion of acceptor states in germanium. These 
have energies of 0.013 ev or less (see Table I), which are much smaller 
than the spin orbit separation X of 0.3 ev. Under these conditions the 
functions Fs and Fe, associated with the split-off band are negligible com- 
pared to Fi, . . , F 4 and the six coupled equations (9.1) may, to a good 
approximation, be replaced by the four equations 




(r) +-Fj(t) = EFjiT) 

KT 


j = 1 , 


4. 


( 10 . 1 ) 


The solutions of (10.1) [or (9.1)] are either fourfold or twofold degenerate, 


At Carnegie Institute of Technology, Pittsburgh, Pennsylvania; for a preliminary 
note see W. Kohn and D. Schechter, Phya. Rev. 99, 1903 (1955). 
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corresponding to the representations Fg (fourfold) and Fs, F? (both two- 
fold) of the tetrahedral double group.®®® 

To obtain approximate solutions of (10.1), we expand the Fj(r) in 
spherical harmonics. Since the effective Hamiltonian of (10.1) is even, 
every solution involves only even or only odd I values. Thus we write 
for even \ for example^ 

2 "lkUr)YiAe,'<>). ( 10 . 2 ) 

2 , , . m 

In the spirit of the variational method we restrict I to values less than 
some U and use a simple trial function for the radial parts /y,j.,n(^). Never- 
theless, the problem would still be quite unmanagable if the form of 
(10.2) could not be determined largely from symmetry considerations. 

As a specific example let us consider the ground state. If we take 
Z < 2, we find that (10.1) has four degenerate solutions, a typical one of 
which has the form 



(10.3) 

where /i, /2, and fz are radial functions. The four rows in the '‘spinors^^ 
represent fi, . . . F4. An approximate eigenfunction and energy is now 
obtained in the usual way by varying the three functions /i(r) to maximize 
the expectation value of the matrix Hamiltonian appearing in (10.1). 
We have used the functional forms 

/i(r) = (10.4) 

and found an ionization energy of 0.0089 ev, with 

ri = 43.3 A, ft — Tz = 33.8 A, 

ai = 1.71 X 10®cm-^a2 = -2.29 X lO^^cm-^aa = 4.97 X lO^^cm-^.®^ 
The inclusion of the neglected functions Fs and Fz by perturbation theory 
increases this energy by less than 1 %. 

The theoretical value of the ionization energy, namely 0.0089 ev, is 
in satisfactory agreement with the measured energies, which range be- 

R. H. Parmenter, Phys. Rev. 100, 573, Table VII (1955); G. Dresselhaus, Phys. 
Rev. 100, 580, Table I (1955). 

W. Kohn and D. Schechter.*® 
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tween 0.0102 and 0.0112 ev (Table I), particularly when one recalls that 
any improvement of the trial function would increase the theoretical value 
slightly — perhaps by 5 to 10%. The size of the remaining discrepancy, 
about 10-20%, is just what might have been expected. For the variation 
of the measured energies from acceptor to acceptor is 10%, which gives 
an estimate of the specific effects arising from the neighborhood of the 
impurity ions. Such effects are of course not included in the effective 
mass theory. 

The excited acceptor states in germanium have not yet been ob- 
served, but the positions of four low-lying levels have been calculated by 
Schechter®® and are listed — together with the ground-state energy — in 
Table VII. These levels are labeled 2p^*^ because their envelope functions 
become hydrogenic 2p functions in the limit where the mass parameters 
(L — M) and N become zero. The superscript (z) merely serves to differen- 
tiate them. It may be remarked that there are 12 of these “2p^^ states in 
all, corresponding to the 12 linearly independent products The nota- 
tion “Is’^ for the ground state has a similar significance. The fourfold 
degeneracy of this state corresponds to the four linearly independent 
functions ^y. 

Table VII. Acceptor Spectra in Effective Mass Approximation 


Energies in units of 0.01 ev 


State Degeneracy Silicon® Germanium 


l9 

4 

3.38 

(4.90) 

0.89 

2pU) 

4 

1.26 

(1-32) 

0.38 

2p'‘> 

4 

0.86 

(0.89) 

0.24 


2 

0.64 

(0.65) 

0.19 

2p«> 

2 

0.43 

(0.41) 



® The first column under ** Silicon corresponds to a negative sign of the mass param- 
eter B s — Af), Eq. (8.6), which is slightly favored by band theory. The second 
column, in parentheses, corresponds to positive B. 

b, Silicon 

It was not an accident that we discussed germanium first, for we en- 
counter many troubles in silicon. 

First of all the mass parameters L, Af, and N in silicon are known only 
very roughly from cyclotron resonance experiments [see Eqs. (8.61)]. 
Even the sign of the combination B = i(L — M) is uncertain, although 
Kane®^ has presented an argument which favors a negative sign slightly. 

Second, the spin orbit splitting X (« 0.035 ev) is now of the same order 
of magnitude as the energies of the acceptor levels («0.05 ev), so that we 
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expect Bloch waves associated with all six valence bands to be involved 
in the acceptor wave functions. Consequently we cannot reduce our six 
coupled equations (9.1) to four as we did for germanium. 

Finally, X has not yet been measured in silicon. The theoretical value 
might be in error by as much as 20%. Fortunately this uncertainty 
affects the calculated energy levels by only about 5%. 

All calculations have been carried out using the values 

B ^ {L - M)/3 = ±1.24ftV2m, 

Both results are listed in Table VII. 

Let us now compare the results with experiment. Beginning with the 
ground state, we see from Table I that, apart from In, the experimental 
ionization energies range from 0.045 to 0.071 ev. This alone shows that 
the specific effects of each type of impurity are important and that an 
effective mass description can have only qualitative significance. Under 
these circumstances the agreement” with the theoretical values of 
0.034 (or 0.049) ev is about what would be expected. The comparison 
between theory and experiment is similar to that for Si donors (Section 
6a). 

In has an ionization energy of 0.16 ev which is about four or five 
times the effective mass value. It is therefore questionable whether it 
should at all be regarded as a shallow” acceptor. For this impurity 
Shulman®^® has proposed a theory which accounts for the observed energy 
qualitatively (within a factor of about 2) in terms of the difference be- 
tween the energies necessary to break an In-Si bond and a Si-Si bond. 

Turning next to the excited acceptor states in silicon, we find a rather 
unclear situation. We recall first that the spacings between the observed 
excited donor levels in silicon were very nearly the same for all donors 
[see Eq. (6.1)]. This showed that the specific effects of the different im- 
purity ions are very small. One would expect the same to be the case for 
the optically excited acceptor states, whose orbits have similar sizes and 
whose envelope functions also have nodes at the position of the impurity 
ions. In our opinion it is not clear from the available experimental data 
whether this is so or not. 

The most accurate and complete data are those of Hrostowski and 
Kaiser,^® shown in Fig. 2. They are consistent with the original work of 
Burstein et aU and of Newman^^® but show considerably more detail. 
The general similarity of the three spectra associated with B, AI, and Ga 
is striking enough. The most disturbing feature is the broad two-peaked 
line for Ga at 500-520 cm“^ As all other lines have much smaller widths, 

R. G. Shulman, private communication. 



SHALLOW IMPURITY STATES IN Si AND Ge 


305 


it appears likely that this structure has, at least in part, a different origin. 
The location of the two B lines at 243 and 278 cin“^ lines has a rather large 
uncertainty of about ±10 cin“^ for instrumental reasons. The *'best'' 
spectrum therefore is that associated with Al. It gives the following 
energy differences between successive excited states: 

Al: 0.37, 0.56, 0.08, 0.21, 0.12 X lO’^ ev. (10.5a) 

For B we have 

B: ^0.5, '-'0.35, 0.12, 0.17, 0.07, 0.06 X lO^^ ev. (10.5b) 

If in the interest of obtaining the best possible agreement with the other 
spectra we arbitrarily discard the peak at 519 cm”^ (0.0643 ev) for Ga, 
we obtain 


Ga: 0.37, 0.58, 0.15, 0.30 X lO-^ ev. (10.5c) 


In view of the foregoing remarks, we can only conclude that the excited 
B and Ga levels may agree with those of Al within about 10%.®^ From 
Table VII we have the following theoretical level spacings: 


. I B < 0 0.40, 0.22, 0.21 X 10“^ ev 
ineory. | g ^ q q 43^ q 24^ 0.24 X 10-^ ev. 


( 10 . 6 ) 


The difference between the lowest two excited levels is in fairly good 
agreement with the experimental results (10.5), The other spacings do 
not agree well. It may be noted, however, that for all three impurities 
the experimental value of the lowest excited state lies 1.3 to 1.4 X 10”^ ev 
below the continuum, which is in very good agreement with the theoret- 
ical estimates of 1.3 X 10“^ ev. In judging the agreement between theory 
and experiment, the uncertainties in the mass parameters and in X, as 
well as the difficulty of obtaining accurate solutions of the six coupled 
equations (9.1), should be remembered. 

We should like to add a remark about the intensities of the absorption 
lines. In atomic hydrogen the transitions from the ground state to the 
2p, 3p, 4p . . . levels have relative intensities of 1.00, 0.19, 0.07, . . . . 
No such preponderance of the lowest transition is observed in the accep- 
tor spectra. Qualitatively, this is in agreement with the fact that in the 
effective mass theory the four lowest lying excited levels all have 
character; that is, all these states become hydrogenic 2p functions when 

A variation of 10% of the absolute energies of the excited states corresponds of 

course to a much longer fractional variation in their energy differences. 
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the mass parameters (L — M) and N are zero,®^ Quantitative theoretical 
calculations of the intensities are, however, very difficult as the required 
matrix elements also involve the ground-state wave functions which are 
only imperfectly known. 

V. Effects of Strains and of Static Electric and Magnetic Fields 

Up to this point we have discussed the structure of donor and acceptor 
states only in a hypothetical perfect crystal, unperturbed by internal or 
external strains, by electric and magnetic fields, or by lattice vibrations. 
In this chapter and the next we shall discuss briefly some of the interest- 
ing effects due to such perturbations. 

11. Strains 

Bulk samples of silicon and germanium can be subjected to elastic 
strains to a maximum of about 3 X 10“®.®^Strainsof theorderof 3 X 10“^ 
are produced fairly easily. There are also internal strains arising from 
dislocation lines whose magnitude is given approximately by 

s « (V^) X 10“® cm. (11.1) 

Here n is the number of dislocation lines per cm^. This gives s « 10“^ for 
the very low dislocation line density of 100 lines/cm^, whereas one has 
s « 10“^ for badly deformed crystals with n = 10® lines/cm^. The effects 
of small strains on the properties of an impurity state associated with a 
single, simple band minimum or maximum are small, of the order of 
magnitude of s itself. However, in silicon and germanium, much larger 
effects should be associated with shear strains because of the degeneracy 
of the conduction band minima and of the valence band maximum. Iso- 
tropic compression still has a negligible effect. So far, the effects of 
anisotropic strain have not been confirmed experimentally. The following 
discussion is intended primarily to indicate possible experiments of 
interest. 

The first experimental spectrum obtained by Burstein and co-workers* for B had a 
fortuitous resemblance to a simple hydrogenic spectrum. Therefore these authors 
labeled their optically excited states in order as 2p, 3p, and 4p. It is now known 
that this labeling has no significance. First of all, Hrostowski^* has discovered an 
additional lower-lying state for B which removes any similarity with the hydrogenic 
spectrum. Also the other acceptor spectra show no such similarity. Finally, from a 
theoretical point of view, no resemblance to the hydrogen spectrum is to be expected. 
In fact, the four lowest strong transitions are all believed to be 2p-like states. 

G. L. Pearson, W. T. Read, and W. L. Feldman, Acta Met. 6, 181 (1957); also W. T. 
Read, private communication. 
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a. Donor States 

We mentioned in Section 4 that the conduction band has several 
equivalent minima in both silicon and germanium, provided the lattice 
is not strained. For definiteness we shall talk, in the subsequent discus- 
sion, about silicon, which possesses six equivalent minima. When the 
crystal is subjected to a specific shear strain of magnitude s, the equiva- 
lence of these minima is destroyed. Some will be lifted and others lowered. 
We can write for the shift of the jth minimum 

A£(ky) = S€,. (11.2) 

Here the e, are proportional to the deformation potential constant Eu 
of Herring and Vogt,®^ and are of the order of magnitude of 6 ev. For 
shear strains, the mean shift of the six minima vanishes: 

6 

J = 0. (11.3) 

3 = 1 

Let us now examine the effect of such a strain on the ground-state 
wave function®® of a donor. In the absence of strains and in the effec- 
tive mass approximation, the ground state is sixfold degenerate. The 
wave functions belonging to different irreducible representations of the 
tetrahedral group are of the form 

6 

^(0 = ^ i = 1, . . . 6 (1L4) 

3 = 1 

[see Eqs. (5.48), (5.45), and (5.49)]. This degeneracy is partly split as a 
result of corrections to the effective mass formalism, but the wave func- 
tions are still given approximately by (11.4). The lowest, nondegenerate 
state is denoted by If the strain s is now produced, first-order per- 
turbation theory shows that the ground-state wave function will be mod- 
ified as follows: 

6 fi 

^,(0) = ^(0, + s ^ (11.5) 

t = 2 3 = 1 

Here the are the energies of the 6 “Is ’Mike states in the absence of 
strain. We see that the change of the wave function depends critically 
on the energy differences It was remarked in Section 5b that 

C. Herring and E. Vogt, Phys. Rev. 101, 944 (1956). 

The importance of such effects has recently been pointed out by P. J. Price, Phys. 

Rev. 104, 1223 (1956). 
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the for i > 0 all vanish at the donor nucleus. This leads one to believe 
that the corresponding will be given more accurately than the energy 
by the effective mass value. Unfortunately, there is no experimental 
evidence concerning the positions of these higher lying Is 'Mike states. 
Strain-induced changes of would throw some interesting light on 
these levels if they could be measured. 

For a rough estimate, let us replace the energy by A, where, 

in accordance with the above remarks, 

A = Eohu JSfeff. maw. (11*0) 


Then (11.5) gives for the approximate fractional change of the wave 
function 


^(0) ~ ® |Af 


(11.7) 


Here c is of the order of magnitude of the €y, that is, about 6 ev. 

For Li in silicon A = 4 X 10~® ev (Tables I and V). Thus a moderate 
strain s — 10”^ produces a 15% change in the wave function — a very 
substantial effect. It must be noted, however, that since all with i > 0 
vanish at the donor nucleus, the first-order change of the wave function 
at this nucleus, given by (11.5) happens to vanish. The second-order 
change is, of course, much smaller. This suggests that one look for changes 
of at other points. The double resonance technique, discussed in 
Section 2b, which yields the values of at the sites of the Si^® nuclei 
within the donor orbit, is particularly well suited for this purpose and 
should yield interesting information. 

For other donors in silicon, A is substantially larger than in the case 
of Li ranging from 9 X 10“® ev in Sb to 4 X 10“^ ev in Bi. The effect of 
strains on the wave function should be proportionately smaller. 

Let us turn next to the change in ionization energy produced by shear 
strains. The first-order change of the ground-state energy is given by 
^ and vanishes, according to Eq. (11.3). This conclusion does 

} 

not depend on the effective mass approximation but follows from general 
symmetry considerations. The second-order change is given by 

( 11 . 8 ) 

where i is again of the order of magnitude of ey. At the same time, the 
strain depresses the lowest point of the conduction band by an amount 
slcmaxl, where is the magnitude of the numerically largest €j. The 
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energy difference Eq (>0) between the lowest point of the conduction 
band and the lowest bound state is, therefore, changed by 

dEo = -s\e^^\+~ (11.9) 

(See Fig. 13.) Equation (11.9) represents the change of the effective 
ionization energy if the temperature is sufficiently low that the occupa- 
tion of all bound states other than the lowest 
bound state and of all conducting states other 
than those near the lowest conduction band 
minimum can be neglected. These conditions 
maybe attainable in favorable cases (e.g., P or 
As in germanium, T » 6°K, s = 3 X 10“^). 

Hall measurements could then give interest- 
ing information about A as well as jemaxl and 
hence provide values of the deformation 
potential constant 

We shall now consider the excited donor 
levels briefly. Each optically active level 
belongs to the vector representation Ti of the 
tetrahedral group and is threefold degenerate 
intheabsenceof strains. A strains = 3 X 10“^ 
may be expected to split this degeneracy by 
roughly sc' « 2 X 10“® ev, and should result 
in an appreciable broadening (or even split- 
ting) of the absorption lines. 

Finally, a remark about internal strains. In ‘‘good" crystals, having 
dislocation densities less than 10® lines/cm^, the effects of internal strains 
(s 10~®) on donor states will be very insignificant. 

6. Acceptor States 

We have seen in Section 10, that in the absence of strains all acceptor 
levels are either twofold or fourfold degenerate. Let us consider first the 
ground state, which is fourfold degenerate according to theory. Under a 
shear strain s, this level will split into two twofold degenerate levels.®® 
The magnitude of the splitting is 

AE = sc. (11.10) 

F. J. Morin, T. H. Geballe, and C. Herring, Phys, Rev. 106, 525 (1956). 

According to Kramers^ theorem (time reversal), at least twofold degeneracy remains 
unless there is an external magnetic field. 



Fig. 13. Effect of a shear 
strain s on the donor ioniza- 
tion energy (schematic). Note 
the linear s dependence of the 
conduction band minima, and 
the quadratic a dependence of 
the ground state level 
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Here e is an effective deformation potential constant, which depends on 
the geometry of the strain and is of the order of 10 ev. The center of 
gravity of the four states remains fixed to first order in s. The wave func- 
tions associated with each of the two levels depend on the nature of the 
strain. 

We see from Eq. (11.1) that internal strains arising from dislocations 
produce a broadening of the lowest acceptor level of the order of 

AE « (10“-® cm) ^/n €. (11.11) 

This is about 3 X 10“^ ev for a typical dislocation density of lOVcm^. 
Such broadening has important consequences for the magnetic properties 
of acceptor states, which are discussed further in Sections 13b and 15c. 



Fig. 14, Effect of a shear strain s on acceptor levels (schematic). Fourfold levels 
are split to first order in s, twofold levels remain unaffected. 

We turn next to the excited levels. Again, the fourfold degenerate 
ones will be split, according to Eq. (11.10), whereas the twofold degen- 
erate ones remain unsplit and fixed to first order in s. 

Let us now see how an external shear strain s = affects the in- 
frared absorption spectrum of acceptor states. At temperatures for which 
kT » se (i.e., T ^ 10°K) both levels which arise from the ground state 
will be occupied equally. Hence all absorption lines should broaden by an 
amount se « 10”^ ev due to this cause. Furthermore, transitions to four- 
fold degenerate final states will be broadened additionally as a result 
of the splitting of the final states. On the other hand, if kT <$C (i.e., 

T 10°K), only the lowest of the two ground-state levels will be 
occupied. Therefore, the transitions to the twofold excited states are 
not broadened by strain; however, the corresponding photon energy is 
increased as a consequence of the depression of the ground-state level. 
Transitions to fourfold excited states will be both broadened, because 
of the splitting of the excited level, and shifted to higher photon ener- 
gies, as a result the movement of the ground state. (See Fig. 14.) Such 
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observations would be interesting for two reasons. They would enable 
one to determine the degeneracy of the various acceptor states experi- 
mentally, and they would give one interesting information about the 
magnitude of the deformation potential constant «. 

Finally, one would expect to obtain interesting effects in strained 
crystals with the use of 'polarized light. For, as we have just seen, only two 
of the four normally degenerate states are occupied when kT <^S€. These 
have preferred orientations'^ relative to the crystal axes. Consequently, 
the intensities of the absorption line depend strongly on the direction of 
polarization of the incident light. 


12, Stark Effect 

The Stark effect, which would be of great interest, has not yet been 
observed in impurity states. However, since it might be detected in sili- 
con, we shall sketch the theory briefly. 

Let 8 be the electric field at the impurity center and let us assume its 
direction is along the z axis. The perturbation Hamiltonian then is 


57 = -^ez, ( 12 . 1 ) 

The wave functions of the bound carrier have the form 


^«(r) 

3 


( 12 . 2 ) 


in the effective mass theory of donor or acceptor states. [See Eqs. (5.48), 
(5.46), and (9.21).] The matrix element of 8V between two such states is 
given by 

^ (12.3) 

3 


Since the effective mass Hamiltonians for both types of states are invari- 
ant under inversion and there are no accidental degeneracies (such as the 
2s, 2p degeneracy in hydrogen) the first-order Stark effect vanishes.®® 
The second-order Stark shift for a nondegenerate state is given by 




‘-1 


|(^( 0 ^ 57 ^( i 0)|2 
^(i) _ £'(i') 


ae ae 


(12.4) 


The full Hamiltonian of the impurity problem has only tetrahedral symmetry and 
is not invariant under inversion. As a result, if the effective mass theory is seriously 
in error, states belonging to the representations Ti, T2, and Tg can have an appreci- 
able first-order Stark effect. This might possibly be of significance for the acceptor 
ground state (Pg) in silicon. 



312 


W. KOHN 


where is an average of the denominators The shifts will be 

of a similar order of magnitude for degenerate states. Let us make a 
rough estimate of h^E, by taking a simple hydrogenic wave function of 
“radius^^ R for This gives 


b^E « 


{eZRY 

AE 


(12.5) 


We see now that impurity states are tempting objects for Stark effect 
measurements because of their large orbits and small energy separations. 

It may be possible to apply about 3000 volts/cm to silicon at low 
temperatures without causing appreciable heating.^® The Stark shift of 
the ground state is, however, only about 3 X 10"“^ ev for such a field, that 
is, is quite negligible. On the other hand, let us consider a highly excited 
state and take as representative values 


R » 10“^® cm 
AF « 5 X 10-3 ev. 


( 12 . 6 ) 


Equation (12.5) then gives 

62F « 5 X 10-3 ev (12.7) 

which is comparable to the unperturbed energies of the excited states. A 
shift of this magnitude could, of course, be detected easily by infrared 
absorption experiments, which have a precision of about 10“® ev (see 
Section lb). 

Stark effect measurements, if feasible, would give us very useful in- 
formation about the excited impurity states. For example, the fourfold 
degenerate acceptor levels should be split into two, whereas the twofold 
degenerate levels should be shifted but not split. 


13. Static Magnetic Susceptibility 

We mentioned in Section 3 that impurity states have been found to 
make significant contributions to the magnetic susceptibilities of the 
samples in which they are located. So far, reliable measurements of the 
contributions of effectively isolated impurity states are not yet available 
but may be expected before very long. It is interesting to see the results 
we would anticipate on the basis of our current models. 


a. Donor States 

In the case of silicon, we know from magnetic resonance measurements 
that the g factor of the donor electron is 1.9985,^^ which is very close to 

I am indebted to Dr. Burstein for a conversation on this question. 

G. Feher, private communication. 
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that of the free electron. This shows that spin-orbit coupling plays a 
negligible role. The corresponding spin paramagnetism is therefore, to 
an excellent approximation, 

2 

Xj. = ^ per donor (13.1) 


where mo is the Bohr magneton. 

In germanium spin resonance experiments have not been successful 
so that we do not have an experimental source for the g factor. The spin 
orbit interaction in germanium may be expected to be about an order of 
magnitude larger than in silicon because of the higher atomic number. 
However, we still expect the free electron paramagnetism (13.1) for all 
practical purposes, provided the ground state has no orbital degeneracy. 

To obtain the diamagnetic susceptibility, one must, as usual, calculate 
the energy of the donor state up the second order in the magnetic field 
H, We have seen that the ground-state wave function is a linear combina- 
tion of several wave packets, each composed of Bloch waves arising from 
the vicinity of one of the conduction band minima. It is easily seen that 
these may be treated independently in the effective mass approximation, 
provided that we eventually take the average of their magnetic energies. 
Since we know from the symmetry of the donor states that the suscep- 
tibility is a scalar, we can also calculate the magnetic energy of any one 
wave packet and average it over all directions of the field. 

The Hamiltonian for one of the packets is 


„„ 1 /h d e .y . I [/h d e . V 

+ (^r~^ /tA] - - ( 13 . 2 ) 

\L dy C / ] KV 


where 


A = i(H X r). 


(13.3) 


To second order in H, the magnetic energy consists of two contributions: 
8Ei, the expectation value of the second-order terms of (13.2); and 3^2, 
the contribution of the first-order terms of (13.2) evaluated by second- 
order perturbation theory. This calculation has been carried out^^ 


using the form 



(13.4) 


as an approximation for the unperturbed ground state. The result is 



e" r j 
6c^ imi rrit 


(a^lmi — b^/nity 
(d^/mi -f- b^/mt) _ 


per donor. 


(13.5) 


SE 2 vanishes for mi = mu 
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With the effective mass values given in Eqs. (4.2) and (4.3) and the char- 
acteristic lengths a and h given in Eq. (6.5), one finds 

Xsi = -2.3 X cgs units/donor 

Xg« = —2.7 X 10”^® cgs units/donor. * 

The fact that the observed ionization energies ^oba are somewhat 
larger than those predicted by the effective mass theory shows that the 
actual orbits are smaller than the theoretical ones. This may be allowed 
for, in a rough way, by multiplying (13.6) by the factor (Ee«.ma«a/Eob«), 
where E^u, masa is 0.029 ev in silicon and 0.0092 ev in germanium. 

6. Acceptor States 

The acceptor ground state is fourfold degenerate (Section 10) in the 
absence of any perturbations. A magnetic field will lift this degeneracy 
completely to first order in H and depress the center of gravity to second 
order. The fixst effect gives rise to a temperature-dependent paramagnet- 
ism, the second to a temperature-independent diamagnetism. 

The paramagnetic susceptibility of germanium has been calculated 
by the writer. The computation is quite lengthy because of the matrix 
nature of the effective mass Hamiltonian. One interesting feature is that 
the susceptibility depends not only on the usual inverse mass parameters 
L, Af, N of Eq. (8.10), but also very critically on another parameter /c, 
introduced by Luttinger^^ in a study of the quantum effects of cyclotron 
resonance. We have estimated that k lies in the range k = (3.5 ± 0.2). 
This gives the paramagnetic susceptibility 

= 

kT (13.7) 

R = 0.35 ± 0.25. 

If K is omitted from the calculation one finds R ^ 35. The small result 
contained in (13.7) is the result of seemingly fortuitous cancellations of 
large terms and may be altered appreciably by a slight improvement of 
the trial wave functions [see Eq. (10.3)]. 

Before applying the result (13.7) to a real sample containing accep- 
tors, the effect of internal strains arising from dislocations must be in- 
cluded. We have seen in Section lib, that a strain s splits the fourfold 
degenerate ground state into two twofold degenerate levels, the separa- 
tion being the order of s X (10 ev). If this splitting is comparable to kT, 
the susceptibility is reduced below the value given by (13.7). The critical 
dislocation density is about 10® dislocation lines/cm^ at T = 1.4®K. 

J. M. Luttinger, Phys, Rev. 102, 1030 (1956) ; this constant k (which occurs nowhere 
else in the present article) has no relation to the dielectric constant. 
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No calculations of the paramagnetic susceptibility of acceptor states 
in silicon are available. The diamagnetic susceptibility of acceptor states 
has not yet been calculated for either silicon or germanium. 


14. Magnetic Field Dependence of the Ionization Energy 

Both the energy of the impurity ground state and of the lowest point 
of the continuum are significantly altered in a strong magnetic field. 
Since the two energy changes are of different magnitudes in general, a 
change of the ionization energy will result. 

For the case of a simple donor state, described by an isotropic mass 
m*, this effect has been calculated by Yafet et alJ^ The strength of the 
magnetic field is conveniently measured by the dimensionless parameter 


ho3l2 


(14.1) 


Here Eo is the ionization in the absence of a magnetic field (the “Rydberg 
constant” of the impurity state) and ho3/2 is the zero point energy of a 
conduction electron of mass m* in the external magnetic field : 


ho) _ eh jj 
2 2m*c 


(14.2) 


The interaction energy of the electron spin with the magnetic field is 
identical in the ground state and the conduction band and can be neg- 
lected. Instead we are concerned here with the effects of the magnetic 
field on the orbital energy of the electron. 

The energy of the edge of the conduction band is simply given by 
(14.2). For weak fields, the energy of the ground state can be calculated 
by perturbation theory, but it must be obtained variationally for stronger 
fields (t 1). Yafet et aL find that the wave function is appreciably 
compressed in all directions for 7 ^ 1, but least of all in the direction 
paralleled to the field. The total energy of the ground state rises, but 
somewhat less rapidly than that of the conduction band edge. The result- 
ing effect is an increase of the ionization energy. The numerical results 
are shown in Fig. 15. 

The magnetic field for which 7 = 1 is given by 

H = 2X lOM— ) oersted. (14.3) 

\Kmo/ 

It is of the order of 10® oersted for the impurity states in germanium and 
is about 10® oersted for those in silicon. 

73 Y. Yafet, R. W. Keyes, and E. N. Adams, J. Phys. Chem. Solids 1, 137 (1956). 
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Fig. 15. Energy difference between the lowest state of the conduction band and 
the impurity ground state as function of the parameter y which characterizes the 
strength of the magnetic field, (After Yafet et aU*) 

The field is only about 1300 oersted in InSb where m^fmo = 0.013 and 
K = 16. Keyes and Sladek^^ have measured the Hall coefficient in this 
material in fields up to 80,000 oersted at 4.2°K and have found a decrease 
of the Hall coefficient by a factor of 50. This large effect can be under- 
stood qualitatively in terms of the change of ionization energy. Quan- 
titative comparison with theory was not possible because of overlap 
effects between neighboring impurities. 

VI. Interaction with Lattice Vibrations 

In this chapter we shall discuss briefly the effects of lattice vibra- 
tions on impurity states. These effects have been studied in considerable 
detail by Lax and Burstein^® in order to account for the observed widths 
of the infrared absorption lines (Section lb). We believe also that the 
failure of the attempts to observe spin resonance of acceptor states may 
arise from the interaction of the electrons with the lattice vibrations. 

15. Absorption Line Widths 

The following discussion is a brief summary of the work of Lax and 
Burstein.’^ Further details may be found in the original papers. 

Let us focus our attention on B acceptor states in silicon for the sake 
of definiteness. First we want to show that it is reasonable to use the 
adiabatic approximation, which is familiar from the theory of molecular 
vibrations. If Eq denotes the ionization energy of the impurity electron, 
a typical electronic frequency is 

Wei = Eo/h « 7 X 10^® sec“^ 

R. W. Keyes and R. J. Sladek, J. Phys, Chem, Solids 1, 143 (1956). 

M. Lax and E. Burstein, Phys. Rev. 100, 592 (1955). 


(15.1) 
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Now the electron interacts appreciably only with acoustic lattice vibra- 
tions whose wavelengths are comparable to the diameter of the orbit. 
Vibrations of much longer wavelength describe a slow rigid displacement 
of the entire impurity state and the interaction with short wavelength 
vibrations is weakened by partial cancellation of the contributions from 
different parts of the impurity orbit. The frequency of the important lat- 
tice vibrations is approximately 


wvib « 1.3 X lO'^sec-i (15.2) 

which is appreciably less than the electronic frequency (15.1). This is the 
condition for the validity of the adiabatic approximation. 

In this approximation, we can write the wave functions describing the 
acceptor ground state and the lattice vibrations in the form 

^„(o) = ^(o)(r)$,(o)(ai, a2, . . .), (15.3) 

Here is the electronic wave function, oi, a 2 , . . . are a set of 

coordinates describing the lattice vibrations and the are the wave 
functions describing various degrees of excitation of the normal modes 
of vibration. A weak dependence of on the a*- has been neglected. 
Similarly, the wave functions corresponding to the ith excited acceptor 
state have the form 


>lrn(0 = ^<')(r)^„(*)(ai, a2, . . .). (15.4) 

It is an essential consequence of the electron-lattice interaction that the 
normal modes are slightly different from 

a. Optical Absorption Line Widths 

The matrix element for an optical transition from to con- 
tains, apart from the electronic dipole matrix element, the scalar product 

(<I>n^®>,<i>n'^'''). (15.5) 

Now if the sets of lattice functions ^n^®^ and were identical (no 
electron lattice interaction), optical transitions could occur only for 
n = w', i.e., between states which do not differ in their vibrational mo- 
tion. The line would then be sharp^^ and given by hv = the 

difference between the electronic energies. Because of the difference be- 
tween the sets 4»n^®^ and however, transitions can take place from a 
given n to a variety of n', that is, transitions in which the vibrational 
state of motion changes. This gives rise to a line width. The shape of 

The width due to the finite lifetime of the excited state is negligible. 
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the absorption line associated with the electronic level i clearly is given by 
n'~>{hv)d{hv) = C ^ 

n,n' 

hv < <hv + d{hv) (15.6) 

in which I is the intensity, C is a normalization constant, and En^^^ and 
En'^*^ are the total energies of the states and respectively. The 
Boltzmann factor describes the occupation probability of the initial state. 

The magnitude of the electron-lattice interaction, which is responsible 
for the difference between the vibrational states associated with and 
can be estimated from the observed lattice mobility of the electrons, 
since electronic wave functions of long wavelengths are involved in both 
cases. In this way Lax and Burstein^^ find that the root-mean-square 
widths of all optical absorption lines are given at T = 0° by 

to a good approximation. Here the symbols have the following meaning: 
M is the mass per unit cell (twice the atomic mass of silicon); v is the 
velocity of sound ; a* is an effective Bohr radius of the impurity ground 
state; is the maximum propagation vector of a lattice vibration; 
and El is an effective deformation potential constant^® of the order 
of 15 ev. Substituting numbers into (15.7), Lax and Burstein obtain 
A(Av) = 3.6 X 10“® ev, in fairly good agreement with the observed widths 
which are 0,5 to 1.0 X 10“® ev. 

The reason why all lines are expected to have nearly the same width is 
as follows. The interaction between electron and lattice vibrations is so 
much weaker in the more extended excited states than in the ground state, 
that the overlap integrals are essentially independent of L 

The experimental line widths are in fact approximately equal for almost 
all lines (see Fig. 2). 

Finally, a word about the temperature dependence of the line width. 
As we mentioned before, the modes of vibration which interact most 
strongly with the electron in its ground state have a wavelength in the 
vicinity of X « ira*. These modes become appreciably excited at a tem- 
perature To given by 

fcTo « ^a>(X). (15.8) 

For a* = 13.5 A this gives To « 94®K. Near this temperature the width 
of the line would be expected to increase appreciably [see the Boltzmann 

J. Bardeen and W. Shockley, Pkys. Rev. 80, 72 (1950). 
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factor in (15.6)]. This prediction agrees with the observed dependence 
of the B line widths upon temperature reported in ref. 75. 

Thus the main features of the widths of the optical absorption line 
appear to be well understood in terms of the electron lattice interaction. 

5. Absence of Paramagnetic Resonance in Acceptor States 

Although careful attempts have been made to find a paramagnetic 
resonance in p-type silicon,®^ similar to the donor state resonance de- 
scribed in Section 2a, the experiments have not been successful. We should 
like to suggest two possible reasons. 

We have seen in Section 15a, that the optical absorption lines of 
impurity states have a width associated with electron-lattice interactions 
of the order of 10“® ev at T = 0°. For quite similar reasons, we would 
expect the paramagnetic absorption lines corresponding to transitions 
between the magnetic acceptor levels to have comparable line widths.’^ 
Since the separation of the magnetic levels is given by = 4 X 10'® ev 
at 9000 mc/sec, we see that such a broadening would destroy the reso- 
nance completely. 

Apart from this effect, there is the broadening due to internal strains. 
Such strains split the ground-state level into two doublets whose mean 
energy difference ZE is given approximately by Eq. (11.11). For a typical 
dislocation line density of 10® lines/cm^, AE » 3 X 10“® ev, just com- 
parable to the magnetic energy at 9000 me. 

In concluding this article, I would like to express my thanks to Dr. 
E. Burstein, Dr. G. Feher, Dr. H. J. Hrostowski, Dr. G. Picus, Dr. 
P. Price, Dr. W. T. Read, and Dr. W. W. Tyler for providing me with 
unpublished results and for numerous helpful comments. 

Appendix: Deep-Lying Impurity Levels 

The shallow impurity levels which are the main subject of this article 
are associated with elements of groups III and V and with Li. Besides 
these shallow levels a great many deeper ones have been identified which 
are introduced by other impurity elements. 

The levels identified clearly in germanium are shown in Fig. 16. Work 
on these levels has been carried out largely by a group at the General 
Electric Research Laboratory.^® These investigators have presented 

This question has not yet been fully analyzed. Preliminary studies indicate a very 
complicated line shape with peaks whose widths may be an order of magnitude 
smaller than the optical width of ev. 

The first deep-lying level in germanium was established by W. C. Dunlap, Jr., 
Phys. Reu. 86, 945 (1952). The review by Burton^ and a recent article by W. W. 
Tyler and H. H. Woodbury, Phys. Rev. 102, 647 (1955) contain references to subse- 
quent work. 
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Fig. 16. Firmly established impurity levels in germanium. All levels are acceptors 
unless marked by D for donor, ka is the distribution coefficient, i.e., the ratio of the 
concentration in the solid phase to that in the liquid phase, (Courtesy W. W. Tyler 
and H. H. Woodbury.) 


strong evidence that the monovalent metals Cu and Au can accept 1, 2, 
or 3 electrons/® whereas the divalent metals Zn and Cd can accept 1 or 2. 
This suggests that the maximum number of acceptable electrons is deter- 
mined, at least in part, by the number required to complete all tetrahe- 
dral bonds. Apart from this feature, the detailed structure of the deep- 
lying states is not yet known. 

In silicon, Au introduces a donor level 0.35 ev above the valence 
band®®’®^ and an acceptor level 0.54 ev below the conduction band.®^ 
Zn introduces an acceptor level 0.31 ev above the valence band.®^ Mn 
produces a donor level 0.53 ev below the conduction band.®® Moreover, 
Fe introduces at 0.45 ev above the valence band and one 0.55 ev below 
the conduction band.®^ Finally, we mention again that one of the group 
III elements, namely In, introduces a fairly deep acceptor level 0.16 ev 
above the valence band. 
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